WELL-POSEDNESS FOR THE EULER-NORDSTROM SYSTEM WITH 
COSMOLOGICAL CONSTANT 



JARED SPECK 



Abstract. In this paper the author considers the motion of a relativistic perfect fluid with self-interaction 
QQ mediated by Nordstrom's scalar theory of gravity. The evolution of the fluid is determined by a quasilinear 

hyperbolic system of PDEs, and a cosmological constant is introduced in order to ensure the existence 
of non-zero constant solutions. Accordingly, the initial value problem for a compact perturbation of an 
^SJ infinitely extended quiet fluid is studied. Although the system is neither symmetric hyperbolic nor strictly 

I i hyperbolic, Christodoulou's constructive results on the existence of energy currents for equations derivable 

from a Lagrangian can be adapted to provide energy currents that can be used in place of the standard 
energy principle available for first-order symmetric hyperbolic systems. After providing such energy currents, 
the author uses them to prove that the Euler-Nordstrom system with a cosmological constant is well-posed 
in a suitable Sobolev space. 

o ! 1. Introduction 

It is well- known that for symmetric hyperbolic systems of PDEs, an energy principle is available that 
implies well-posedness (local existence, uniqueness, and continuous dependence on initial data) for initial 
data belonging to an appropriate Sobolev space. Consult [5], [ID]) [13]) [H]> E3- or [3T] for the definition 
of a symmetric hyperbolic system and a detailed proof of local existence in this case. A full proof of 
well-posedness is difficult to locate in the literature, but Kato [TB] supplies one using a very general setup 
that applies to symmetric hyperbolic systems in a Banach space. Additionally, for strictly hyperbolic (not 
necessarily symmetric) systems, well-posedness follows from the availability of a generalization of the energy 
principle for symmetric hyperbolic systems. For strictly hyperbolic systems, there are a variety of methods 
due to Petrovskii, Leray, Garding, and Calderon for generating energy estimates; consult [9] or |21j for details 
C^l on these methods. 

^vq We consider here the Cauchy problem for the Lorentz covariant Euler-Nordstrom (EN) system, which is a 

scalar caricature of the general covariant Euler- Einstein system describing a gravitationally self-interacting 
fluid. The EN system is a quasilinear hyperbolic system of PDEs that is not manifestly symmetric hyperbolic. 
Moreover, because of the repeated factors in the expression for Q(x; •) in equation ( 5.1.7[ ) below, and because 
the sheets of the characteristic subset of the cotangent space at x intersect (see Fig. [IJ, it is not strictly 

!^ hyperbolic. Therefore, well-posedness for the EN system does not follow from either of these two well-known 

frameworks. 

Fortunately, alternate techniques recently developed by Christodoulou [B], and which are applied to the 
study of relativistic fluids in Minkowski spacetime in particular in [7], offer a viable approach to studying 
the Cauchy problem for the EN system. The central advantage afforded by Christodolou's techniques, 
which provide energy currents for equations derivable from a Lagrangian, is that they bypass the physically 
artificial requirement of symmetry in the equations: even though the EN system is not manifestly symmetric, 
its energy currents allow for precisely the same energy estimates to be made as in the theory of symmetric 
hyperbolic systems. Once one has these estimates, the proof of well-posedness for the EN system mirrors the 
well known proof for symmetric hyperbolic systems. Our main goal is to use the method of energy currents 
to prove the following theorem (stated loosely here), which is divided into parts and stated rigorously in 
Section 

Main Theorem (Well-Posedness). Let N > 3 be an integer. Assume that the initial data 
V for the EN system are an H N perturbation of a constant background solution V. Then these 
data launch a unique solution V possessing the regularity property V — V € C^([0,T] x M 3 ) n 
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C°([0,T],H N ) n C l ([Q,T],H N - x ). Furthermore, the map from the initial perturbation V - V to 
V — V is a continuous map from an open subset of H N into C°([0, T], H N ). 

While Christodoulou's methods are not the only techniques available for proving the well-posedness of 
the EN system, they are powerful and natural in the sense that they exploit the inherent geometry of 
the equations. In contrast, one may proceed by seeking a change of state-space variables that renders the 
system symmetric hyperbolic. For example, Makino applies this symmetrizing technique to the Euler-Poisson 
equations in |23] , and Makino and Ukai apply it to the relativistic Euler equations without gravitational 
interaction in [24 and |25j . Further discussion of applications of symmetrization discussed in the literature 
can be found in sections 3.1 and 3.2| a Yet the symmetrizing method is not without disadvantages: one must 



solve a formally over-determined system of equations to find the symmetrizing variables^] and the resulting 
state-space variables, if they exist, may place un-physical and/or mathematically unappealing restrictions 
on the function spaces with which one would like to work. However, it should be noted that Makino's 
symmetrization is currently capable of dealing with a restricted class of compactly supported data, while the 



techniques applied here cannot yet handle such data due to singularities in the energy current (5.5.1 1 when 
the proper energy density p of the fluid vanishes. 

2. Remarks on the Notation 

We introduce here some notation that is used throughout this article, some of which is non-standard. We 
assume that the reader is familiar with standard notation for the L p spaces and the Sobolev spaces H k . 
Unless otherwise stated, the symbols L p and H k refer to L' P (M. 3 ) and H k (R 3 ) respectively. 

2.1. Notation and assumptions regarding spacetime. In the Euler-Poisson system with cosmological 
constant introduced below, we use tgRto denote the time variable and s e K 3 to denote the space variable. 
In the Euler-Einstein and EN systems (which we also equip with a cosmological constant below) , we assume 
that spacetime is a 4-dimensional, time-orientable Lorentzian manifold M. and use the notation 

(2.1.1) x = (x ,x\x 2 ,x 3 ) 

to denote spacetime points. For the EN system with cosmological constant, we assume the existence of 
a global system of rectangular coordinates (an inertial frame), and for this preferred time-space splitting, 



we identify t — x° with time and s — (x 1 , x 2 , a; 3 ) with space and use the notation (2.1.1) to denote the 
components of x relative to this fixed coordinate system. 

2.2. Notation regarding differential operators. If F is a scalar or finite-dimensional array-valued func- 
tion on IR 1+3 , then DF denotes the array consisting of all first-order spacetime partial derivatives (including 
the partial derivative with respect to time) of every component of F, while V^F denotes the array of con- 
sisting of all a th order spatial partial derivatives of every component of F; this should not be confused with 
V, which represents covariant differentiation. 

2.3. Index conventions. We adopt Einstein's notation that repeated Latin indices are summed from 1 to 3, 
while repeated Greek indices are summed from to 3. Indices are raised and lowered using a spacetime metric, 
which varies according to context. 

2.4. Notation regarding norms and function spaces. If V is a constant array, we use the notation 

\\ F \\lI (A) = \\ F ~ V\\l*(A), 

and we denote the set of all Lebesgue measurable functions F such that 

ll-^IU p - (^i) < 00 by L P ^ / .(A). Unless we indicate otherwise, we assume that A = M 3 when the set A is not 
explicitly written. 

If F is a map from the interval [Ti,T2] into the normed function space X, we use the notation 

III F |IU,pi,T 2 ] = sup \\F(t)\\ X - 
tepi,T 2 ] 



a The references given are far from exhaustive; we merely wish to provide the reader with some examples of the application of 
well-known techniques. 

^Consult chapter 3 of |10| for a discussion of symmetrization. 
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We often abbreviate ||| F \ \\x,t in place of ||| F |||x,[o,t] • 

We also use the notation C fe ([0,T], X) to denote the set of fc-times continuously differentiable maps from 
(0, T) into X that, together with their derivatives up to order k, extend continuously to [0,T]. 

If A C M. d (d frequently equals 3,4, or 10 in this article) and A is open, then C£(A) denotes the set 
fc— times continuously differentiable functions (either scalar or array-valued, depending on context) on A 
with bounded derivatives up to order k that extend continuously to the closure of A. The norm of a function 
F € C%{A) is defined by 

\F\ k>A ^ ]T sup\d s F(z)\, 
\s\<k zeA 

where ds represents differentiation with respect to the arguments z of F (which may be spacetime variables 
or state-space variables, depending on the context). 

2.5. Notation regarding operators. If X and Y are normed function spaces, then C(X, Y) denotes the 
set of bounded linear maps from X to Y. If U G C(X, Y), then we denote its operator norm by ||W||x,y- If 
X = Y, we write C(X) instead of C(X, Y) and \\U\\x instead of ||W||x,x- If W(t, t') is an operator-valued map 
from the triangle 

At = {0 < t' < t <T} into C(X), then we adopt the notation 

III U |||x,A T = f sup \\U(t,t')\\ x . 

(t,t')eA T 

2.6. Notation regarding constants. We use the symbol C to denote a generic constant in the estimates 
below which is free to vary from line to line. If the constant depends on quantities such as real numbers N' , 
subsets A of K d , functions F of the state-space variables, etc., that are peripheral to the argument at hand, 
we sometimes indicate this dependence by writing C(N',A, F), etc. We frequently omit the dependence 
of C on functions of the state-space variables below in order to conserve space, but we explicitly show the 
dependence when it is (in our judgment) illuminating. Occasionally, we shall use additional symbols such as 
Cq 2 ,L,K, etc., to denote constants that play a distinguished role in the discussion below. 

3. The EN and EN k Models in Context 

The EN system is an intermediate model in between the Galilean covariant Euler-Poisson (EP) and the 
general covariant Euler-Einstein (EE) systems for self-gravitating classical fluids. Although it is the most 
fundamental of these models for self-gravitating Eulerian fluids, the EE system presents numerous technical 
difficulties that make a detailed analysis of the system's evolution, through either numerical or analytical 
methods, extremely difficult. For example, in General Relativity there is a coordinate gauge freedom due to 
the diffcomorphism covariance of the equations, and furthermore, there is no known law of local conservation 
of gravitational energy. Our main motivations for studying the EN system are to bridge the gap between 
the EP and the EE systems and to provide a special relativistic primer for studying the EE system. 

Since it is based on Nordstrom's theory of gravity, it should be stressed that the EN system is physically 
wrong. However, since both the EN and the EE systems are relativistic generalizations of the EP system, we 
expect, at least in some limiting cases, that there are some qualitative similarities between solutions to the 
three systems. Furthermore, in |32] . Shapiro and Teukolsky discuss numerical simulations of the EN system 
in the spherically symmetric case; they expect that the numerical schemes developed in their paper can be 
adapted to allow for the calculation of accurate wave forms in the EE model. 

Before discussing the EN system in detail, we briefly recall the EP and EE systems, endowing both with 
a cosmological cons tan^ denoted by n 2 . We also briefly discuss some local existence proofs for these systems 
in the case k = 0, emphasizing their dependence on the symmetric hyperbolic setup or the method of Leray 
(strict) hyperbolicity. 

We introduce a positive cosmological constant out of mathematical necessity: the EN system fails to 
have non-zero constant solutions without it. Our reasoning is similar to the reasoning that led Einstein 
to introduce the cosmological constant into General Relativity; he sought a static universe, and General 
Relativity without a cosmological constant features only Minkowski space as a static homogeneous solution 



c We deviate from Einstein's notation; he denoted the cosmological constant by A. 
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(see [H]). We emphasize the presence of the cosmological constant n 2 in the models by referring to them as 
the EP re , EE K , and EN K systems; note that EP=EPo, and similarly for the other two models. 

3.1. The Euler-Poisson system with cosmological constant (EP K ). In units with Newton's universal 
gravitational constant equal to 1, the equations governing the dynamics in this case are 

d t r/ + v k d k i] = 



(3.1.1) 
(3.1.2) 
(3.1.3) 
where 
(3.1.4) 
and 
(3.1.5) 



d tP + d k (pv k ) = 

p {d tVj + v k d k vi) + d jP + pdrf = (j = 1, 2, 3), 



Acj) - n 2 cf) = Anp 



p = y(p,Ti)- 



The unknowns in (3.1.11 



(3.1.4) are the cosmological Newtonian gravitational scalar potential <f>(t,s), 
and the state-space variables mass density p(i,s), velocity v(i, s) = (vi, V2,vs), pressure p(t,s), and entropy 
densitjj^] rj(t, s). We remark that in the EP K system, <p is not a state-space variable because it is uniquely 
determined by p under the assumption of appropriate decay conditions on <f> and p at infinity. The equation 
that specifies p as a function O 3 of p and r\ is known as the equation of state. 

This system of equations^] is discussed in |19j . in which, under an isothermal equation of state (p = c^p, 
where the constant c s denotes the speed of sound) , Kiessling derives the Jeans dispersion relation that arises 
from linearizing (3.1.21 - ( |3.1.4[ ) about a static state in which the background mass density p is non-zero, 
followed by taking the limit k — > 0. 

In 23j, Makino studies the Cauchy problem for the EPo system^ with "tame" compactly supported initial 
data belonging to an appropriate Sobolev space. He studies adiabatic equations of state (p — Kp" 1 , where 
if is a positive constant) under the mathematical assumption 1 < 7 < 3, and after finding symmetrizing 
variables, he proves local existence using the symmetric hyperbolic setup. 

Remark 3.1.1. Let us now make a few remarks about the "tame" data. Vanishing mass densities typically 
produce singularities in the expression for the energy, but Makino's choice of symmetrizing variables, which 
works for the class of adiabatic equations of state described in the previous paragraph, allows him to handle 
a class of compactly supported data. The "tame" data are constrained by the requirement that p s must 
belong to an appropriate Sobolev space, where S is a positive constant depending on 7. To the author's 
knowledge, a fully satisfactory treatment (i.e., without unphysical mathematical restrictions on the data) of 
the evolution of compactly supported data in the EPo system remains an open problem. 

3.2. The Euler-Einstein system with cosmological constant (EE K ). We work in units with Newton's 
universal gravitational constant and the speed of light both equal to 1. Given T, the energy- momentum tensor 
of the contemplated matter model, the gravitational spacetime with cosmological constant is determined by 
the Einstein field equations, 



(3.2.1) 



87rT, 



ft V 



(0 < (i,u<3), 



where G is the Einstein tensor of the spacetime metric g. As a consequence of (3.2.1 ), T has to satisfy the 
admissibility condition 



(3.2.2) 



V M T^ = 



(0 < v < 3), 



where the V denotes the covariant derivative induced by the spacetime metric g. Equation (3.2.2) follows 
from the twice contracted Bianchi identity, which implies that 



(3.2.3) 



0. 



We are influenced by Boltzmann's notation in denoting the entropy density by rj. 
°Kiessling omits equation ( | 3- 1 . 1 [ ) from the system of equations he studies. See Section 
^Equation | |3.1.1| is also omitted from Makino's paper. 



3.2 



for further discussion of this truncation. 
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together with 
(3.2.4) 







(0< A,/x,i/<3), 



which follows from the fact that V is the Levi-Civita connection on spacetime. 

We now briefly introduce the notion of a relativistic perfect fluid. Readers may consult pQ or [H] for more 
background. For a perfect fluid model, the components of the energy-momentum tensor of matter read 

(3.2.5) = (p + p)u^u v +pg^ u . 

Here the scalar p > is the proper energy density, the scalar p > is the pressure, and the vector u is the 
^-velocity, a future-directed timelike vectorfield which is subject to the normalization condition 

(3.2.6) g^u" = -1. 

We also introduce the additional thermodynamic scalar variables n > 0, the proper number density, and 
rj > 0, the proper entropy density, and the following continuity equation: 



(3.2.7) 



V^rm") = 0. 



When g is given and T is defined by (3.2.5 1, equations (3.2.2 1 and (3.2.71 together form the Euler equations 



for a gcncral-rclativistic perfect fluid. In general, when both g and T are unknowns, (3.2.1 I, its consequence 



(3.2.21, and (3.2.5) - (3.2.71 form the EE K system for u, p,p>,n,r\, and g (up to closure, for instance by 
providing two equations that relate p,p,n, and r/). As in the EP K system, the under-determined system 
consisting of (3.2.1), (3.2.21, and (3.2.51 - ( |3.2.7 ) may be closed (up to a choice of coordinate gauge) by 
providing further relationships between the state-space variables. An example of a simple closure often 
discussed in the mathematical (consult e.g. [T], |24j, [25]) and astrophysical (consult e.g. [30]) literature is to 



assume that p is a function of n alone, in which case equation (3.2.71 is an automatic consequence of (3.2.2 1, 
(3.2.51 and the thermodynamic relation ( |3.3.4| below. Equivalently, one may specify p as a function of p 
alone; such fluids are called barotropic. If the fluid is barotropic, the variable rj becomes passive in the sense 
that it satisfies the equation u^V ^rj = 0, but does not otherwise enter into the dynamics; the remaining 
state-space variables (which we may take to be u,g,p) decouple from rj. 

Local existence for a closed relativistic fluid system has been discussed by several authors under various 
assumptions. For example, in 5 , Choquct-Bruhat showed that the EEo system with pressure-free dust 
sourcesp] forms a well-posed Leray-hyperbolic system, and in |30j . Rendall adapted Makino's symmetrization 



(as discussed in Section 3.1 ) of the EP system to handle a subclass of compactly supported initial data for 
the EEo system with perfect fluid sources under an adiabatic equation of state with 7 > 1 . Similar results 
are also proved in [2J, in which Brauer and Karp write the equations as a symmetric hyperbolic system in 
harmonic coordinates. 



3.3. The Euler-Nordstrom system with cosmological constant. We bas e o ur d iscussion here on 
Calogero's derivation of the Nordstrom- Vlasov system^ [3]. Consult sections 



2.1 



and 



2.3 



for some remarks 



on our assumptions concerning spacetime and our use of index notation. As in the EE K model, we work in 
units with the speed of light and Newton's universal gravitational constant both equal to 1. 



Like the EE K system, the EN K system subsumes equations (3.2.21, (3.2.51, (3.2.6), and (3.2.7), where 



P > 0,p > 0, n > Q, t] > 0, and u are defined as in the EE K system. In contrast to the EE K model, we do not 



assume Einstein's field equations (3.2.11; instead we turn to Nordstrom's theory of gravity. We postulate 



that in our global rectangular coordinate system, the conformally flat metric is given by 
(3.3.1) 9„u = e 2 % u , 

where is the Nordstrom scalar potential, and g = diag(— 1, 1, 1, 1) are the components of the Minkowski 
metric in the rectangular coordinate system. 

Nordstrom's theory of gravity [28] belongs to the class of theories known as scalar metric theories of 
gravity. For theories in this class, gravitational forces are mediated by a scalar field (or "potential") (j> that 
affects the spacetime metric. Furthermore, it is assumed that the effect of 4> is to modify the otherwise 



SThe energy-momentum tensor for pressure-free dust has components T M " = pu^u" . 

^Each of the three Eulerian fluid models discussed in this article has a kinetic theory counterpart. Collectively known as the 
Vlasov models, these diffeo-integral systems describe a particle density function / on physical space X momentum space that 
evolves due to gravitational self-interaction. In particular, the ENo system is the Eulerian counterpart of the previously studied 
Nordstrom- Vlasov (NV) system (which does not feature a cosmological constant). See e.g., [3] or 
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flat metric by a scaling factor that depends on <fi. Therefore, the physical metric in such a theory is given 
by g^v = X 2 (4>) 9 , where g is the Minkowski metric. A metric of this form is said to be conformally flat. 
Strictly speaking, the scalar theory of gravity we study in this paper is not identical to the one published 
by Nordstrom in |28j . In his paper, Nordstrom makes the choice x('P) — 4>i while in our paper, we make the 
choice = e^, a theory that appears as a homework exercise in the well-known text "Gravitation" by 
Misner, Thorne, and Wheeler [35]. See [3] or [TT] concerning the significance of the choice x{4>) — which 
has the property of scale invariance of the gravitational interaction. Also consult |29j for a discussion of 
scalar theories of gravity, including the two mentioned here. 

Following Nordstrom's lead [28 , we also introduce the auxiliary energy-momentum tensor T aux with 
components 



(3.3.2) 

and postulate that is a solution to 



(3.3.3) 



□0 - « 2 = -9auTZ = -e 40 (3p - P ). 



Note that □</> = — 9 t 2 </> + A(f> is the wave operator on flat spacetime applied to <f>. The virtue of the postulate 



(3.3.31 is that it provides us with continuity equations for an energy-momentum tensor in Minkowski space 



which we label 6 and discuss below; see equations (4.1.81 and (4.1.91 



As in the EP K and EE K models, we may close the EN K system by supplying relationships between the 
state-space variables. The basic postulates we adopt are as follows (see e.g. [2]): 

1) p > is a function of n > and rj > 0. 

2) p > is defined by 

dp 



(3.3.4) 



def 

p = n 



On 



where the notation |. indicates partial differentiation with • held constant. 
3) A perfect fluid satisfies 



(3.3.5) 



(3.3.6) 



(3.3.7) 



dp 

dn 



>0, 



dp 
dn 



> 0. 



dp 

drj 



> with " 



iff 7/ = 0. 



As a consequence, we have that a, the speed of sound in the fluid, is always real: 



2 45_ f 



dp 
d~ P 



dp/dn\ n 
dp/dn\ v 



> 0. 



4) We also demand that the speed of sound is positive and less than the speed of light whenever 
n, rj > 0: 

< a < 1. 



Postulates 1-3 express the laws of thermodynamics and fundamental thermodynamic assumptions, while 
as discussed in detail in Section [5] postulate 4 ensures that vectors that are timelike with respect to the 
sound cone are necessarily timclike with respect to the light cone. 

Remark 3.3.1. We note that the assumptions p > 0,p > together imply that the energy momentum 



tensor (3.2.5) satisfies both the weak energy condition (T^X^X 1 ' > holds whenever X is future-directed 
and timelike) and the strong energy condition ([T^ u — l/2g a ^ 3 T a pg lll ,\X il X 1 ' > holds whenever X is future- 
directed and timelike). Furthermore, if we assume that the equation of state is such that p = when p = 0, 
then ( |3.3.7 ) guarantees that p < p. It is then easy to check that < p < p implies the dominant energy 
condition (—T^X 11 is future-directed and causal whenever X is future-directed and causal). 



Remark 3.3.2. By (3.3.51, we can solve for a and p as functions of p and r\ 

(3.3.8) <t = Hv,p) 

(3.3.9) P = »(»?, P). 
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Remark 3.3.3. We will make use of the following identity implied by (3.3.61, (3.3.81, and (3.3.91: 



(3.3.10) 



dp 



07, P) 



As a typical example, we mention a polytropic equation of state, that is, an equation of state of the form 
(see e.g. [14"] ) 

(3.3.11) p = n' A{V) ^ 



7 



1 



where 1 < 7 < 2, and A is a positive, increasing function of rj. In this case p = An 1 , dp/dp\ rj is increasing 
in p, and the speed of sound a is bounded from above by y/j — 1. 

Remark 3.3.4. We note here a curious discrepancy that arises when, for the polytropic equation of state 
under the isentropic condition 77 = 770, we consider the Newtonian limit, that is, the limit as the speed 

A c {rj)n\ 



of light c goes to 00. In dimensional units, (3.3.111 becomes p 



m c n 



7-1 



and p 



where mo is the mass per fluid element, and A c (rj) is a positive, increasing function of r\ indexed by the 



parameter c. The speed of sound squared is given by a 



o def n 

z — <r 



7 c J A c ())o)«-'- 1 
c 2 m + (7/7-I) A a (jjo)r 



Assuming that 



lim c _>oo A c (r)o) = j4oo(?7o) exists, we may consider the Newtonian limit c — > 00 of a 2 and p, obtaining in the 
limit that a 2 = "frriQ Aqo^o)?! 7-1 and p = A^-q^iri 1 , Newtonian formulas that make mathematical sense 
and have physical interpretations for 1 < 7 < 00. In the Newtonian case, 7 = 1 corresponds to isothermal 
conditions, while 7 — > 00 yields the rigid body dynamics. However, for finite values of c, not all values of the 
parameter 7 make mathematical or physical sense: there is a mathematical singularity in the formula for p 
at 7 = 1 . This is physically reasonable since isothermal conditions require the instantaneous transfer of heat 
energy. Thus, for finite c, the polytropic equations of state do not allow for the case corresponding to the 
instantaneous transfer of heat energy over finite distances, a feature which we find desirable in a relativistic 
model. Additionally, we have that linin^oo a 2 = c 2 (j — 1), so that for 7 > 2, there is a 7— dependent critical 
threshold for the number density above which the speed of sound exceeds the speed of light. Since larger 
values of 7 correspond to "increasing rigidity" of the fluid, and the concept of rigidity violates the spirit of 
the framework of relativity, we are not surprised to discover that large values of 7 may lead to supcrluminal 
sound speeds. However, we find ourselves at the moment unable to attach a physical interpretation to the 
fact that the mathematical borderline case is 7 — 2. 



We summarize this section by stating that equations (3.2.2), (3.2.51, (3.2.61, (3.2.71, (3.3.11, (3.3.21, 



(3.3.31, (3.3.41, and (3.3.9) constitute the EN K system. 



4. Reformulation of the EN K System, the Linearized EN K System, and the Equations of 

Variation 

Because it is mathematically advantageous, in this section we reformulate the EN K system as a fixed- 
background theory in flat Minkowski space. This is a mathematical reformulation only; the "physical" metric 



in the EN K system is g from (3.3.1) rather than the Minkowski metric g. We also discuss the linearization 



of the EN K system and the related equations of variation, systems that are central to the well-posedness 
arguments. 

4.1. Reformulating the EN K system. For the remainder of this article, indices are raised and lowered 
with the Minkowski metric, so for example, d x 4> = g^ x d^(f>. To begin, we use the form of the metric (3.3.1 1 



to compute that in our fixed rectangular coordinate system (see Section 2.1), the continuity equation (3.2.21 



for the energy-momentum tensor (3.2.5) is given by 

= V^T^ = d^ v + 6T^0 - e-^g a0 T a ^d v (j) 
(4.1.1) =d^ + 6T^d^-e~^g aP T«£d»cb {y = 0, 1, 2, 3), 



where is given by (3.3.21. For this calculation we made use of the explicit form of the Christoffel symbols 



in our rectangular coordinate system: 



(4.1.2) 
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Under the postulate (3.3.31 for 0, (4.1.11 can be rewritten as 



(4.1.3) 



Equation (4.1.31 now illustrates the divergence-free energy-momentum tensor O mentioned in Section 3.3 



Its components Ml/ consist of the terms from (4.1.3) that are inside the parentheses; we are thus afforded 
with local conservation laws in Minkowski space. 



To simplify the notation, we make the change of state-space variables (recalling equation (3.3.91 for the 
definition of the function 3V) 



R = e^p : 
P = e^p 



(i/ = 0,1,2,3) 



(4.1.4) 
(4.1.5) 
(4.1.6) 

throughout the EN K system, noting that U is subject to the constraint 
(4.1.7) U° = (l + U k U k ) 1/2 . 

Following the above substitutions, has components 

(4.1.8) 



and (4.1.3) becomes 
(4.1.9) 



= 



(^ = 0,1,2,3). 



We perform the same changes of variables in the equation (3.2.71 and expand the covariant differentiation 



in terms of coordinate derivatives and the Christoffel symbols (4.1.21, arriving at the equation 

Mu^) = 0. 



(4.1.10) 



For our purposes below, we take as our equations the projections of (4.1.9 ) onto the orthogonal complement 



of U and in the direction of U. In this formulation, the mathematical form of the EN K system is that of the 
relativistic Eulcr equations in Mikowski space without gravitational interaction (as presented in |7J), with 
inhomogeneous terms involving D<j>, and supplemented by the linear Klein-Gordon equation ( 3.3.3[ ) for <fi. 
Thus, we introduce II, the projection onto the orthogonal complement of U, given by 



(4.1.11) 



n"" = WU" + g^. 



Considering first the projection of (4.1.9) in the direction of U, we remark that one may use (3.3.41 and 
(4.1.10) to conclude that for C 1 solutions, U v d jJl @ tl1 ' — is equivalent to 

(4.1.12) U^8^ = 0, 

which implies that the entropy density 77 is constant along the integral curves of U. 



The projection of (4.1.9) onto the orthogonal complement of U gives the 4 equations (only 3 of which are 
independent) 



(4.1.13) 



(R + P)U"d^U" + IP^P = -(□<£ - k 2 0)11^0 (y = 0, 1, 2, 3). 



By ( |3.3.9 ), (4.1.51 and (4.1.6), we may solve for R as a function 91 of rj,P and <f> : 
(4.1.14) R = m(r),P,(l>) = e 4 ^(7 7 ,e- 4 ^P). 

We also the nameless quantity Q and make use of (3.3.41, (3.3.61, ( |3.3.8 ), (3.3.9), (3.3.10), (4.1.51, and 
(4.1.6 1 to express it as a function Q of r],P and <j> : 

* dP 



(4.1.15) 



dn 



v,<t> 



dP 



v,4> 



dp 

dn 



e 4 *S 2 (r,,p)(p + p) 



= S 2 ( V ,e-^P){X(r 1 ,P,0) + P}. 
We also we use the chain rule together with (4.1.10), (4.1.121, and ( 4.1.15[ ) to derive 
(4.1.16) U^P + Qd^U" = (4P - 3Q)U»d^, 
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which we may use in place of (4.1.101. 

Deleting the redundant equation from (4.1.13), using (4.1.71 to derive the relation 

(4.1.17) 



and rewriting (3.3.3) as an equivalent first order system, the working form of the EN K system that we adopt 
is 



(4.1.18) lP i d„r) = 

(4.1.19) U^P + Q^d Q U k + QO k U k = (4P - 3Q)U^ 

(4.1.20) (R + PWdnU* + IP^P = (3P - P)IP J '^ (j = 1, 2, 3) 

(4.1.21) -d a ip + d 3 ipj = k 2 ^ + R - 3P 

(4.1.22) do^- 0^0 = (j = 1,2,3) 

(4.1.23) 9o0 = ^o- 

Here, U°,R, and Q are expressed in terms of the unknowns through the relations 

(4.1.24) U° = (1 + U k U k ) 1/2 

(4.1.25) Q = 0(77,P» 

(4.1.26) fl = lH(» 7 ,P,^) ) 



where the function Q is defined in (4.1.151, and the function is defined in (14.1.14 1. In our rewriting of 



(3.3.3 1 as a first order system, we treat ip„ == d v <f) as separate unknowns for v = 0, 1,2,3. 
To simplify the notation, we collect the unknowns V together into an arrajj^] given by 

(4.1.27) V d = f (77, P, f/\ f/ 2 , (7 3 , 0,^, ^1, V>2, V> 3 ), 
and we refer to the first five components of V as 

(4.1.28) W = (r,,P,U\U 2 ,U*). 

4.2. Linearization and the Equations of Variation (EOV). The standard techniques for proving well- 
posedness require the linearization of the EN K system around a known background solution, which we refer 
to as a "bgs." Each bgs V : A4 — > K 10 we consider is of the form V — (77, P, ■ ■ ■ , ^2,^3)- The resulting system 
is known as the equations of variation (EOV). Thus, given such a V and inhomogeneous terms /,<?,••• , l^ 4 \ 
we define the EOV by 



(4.2.1) 




= .f 


(4.2.2) 


U^d^P + Q^ ) doLf k + Qd k U k 


= 9 


(4.2.3) 


{R+P)U i "d tl U j +U^d M P 


= h® 


(4.2.4) 


-d tpo + d 3 ipj 


= l<°) 


(4.2.5) 


d ipj - djipo 


= 1<D 


(4.2.6) 




= 1^, 



= 1,2,3) 
= 1,2,3) 



'Although every array appearing in this article is a q X 1 column vector, we write them as if they were row vectors to save space. 
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where 




(4.2.7) 


JjO 4£ f 


(4.2.8) 




(4.2.9) 


Q = 


(4.2.10) 


R = 



a/2 



Here, the function is defined in (4.1.15), and the function d\ is defined in (4.1.141. The unknowns 
are the components of V == (fj,P,--- ,^2,^3), and we label the first five components of V by W = f 
(f,,P,U\U 2 ,U 3 ). 

The EOV play multiple roles in this article. Except when discussing the space of variations V as an 
abstract vector space isomorphic to M 10 , we use the symbol V to represent a quantity that solves the EOV. 
The quantity represented by V, the bgs V, and the inhomogeneous terms will vary from application to 
application, but we will always be clear about their definitions in the relevant sections. 

In the case that we are discussing the linearization of the EN K system around a bgs V, the inhomogeneous 
terms take the form 



(4.2.11) 


/ 


= $(V) = 


(4.2.12) 


9 


= (9(V) = (4P - 3Q)[>Vv 


(4.2.13) 




= & j) (V) = (3P - fi)fi«^ 


(4.2.14) 


1(0) 


= £(°>(V) d =V^+i?-3P 


(4.2.15) 


lU) 


= £«)(v)= f (j = 1,2,3) 


(4.2.16) 


;(4) 


= £( 4 )(V) d = f ^ , 



(j = 1,2,3) 



where 5,0,- •• , £ (4) are functions of V. 



It is quite important that the coordinate derivatives of solutions to (4.2.11 - (4.2.10) also satisfy (4.2.11 



(4.2.101 with different inhomogeneous terms. This may be seen by differentiating the equations and 



relegating all but the principal terms to the right-hand side. Similarly, the difference of two solutions to 



(4.2.1 1 - (4.2.101 also satisfies (4.2.1 1 - (4.2.101. Thus, the "•" is a suggestive placeholder that will frequently 



represent "derivative" or "difference" depending on the application. 



Notation. In reference to the inhomogeneous terms on the right-hand side of (4.2.11 ) - (4.2.16 1, we often use 
vector notation including but not limited to 



(4.2.17) 
(4.2.18) 



I =(1^,1^,1^,1^,1^). 



When it is convenient, we will use different vector notation to refer to the inhomogeneous terms, but we 
always use the notation /, g, ■ ■ ■ ,1 to refer to the inhomogeneous terms in scalar form; our use of notation 
for the inhomogeneous terms will always be made clear in the relevant sections. 



Terminology: If V is a solution to the system (4.2.11 - (4.2.101, we say that V is a solution to the EOV 



defined by the bgs V with inhomogeneous terms (b, 1). 

Notation. We will often find it advantageous to abbreviate the "upper half of the various systems in this 



article using matrix notation. For example, we sometimes write (4.2.11 - (4.2.31 as 



(4.2.19) 



A"(V)0 M W = b, 
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where each A M (V) is a 5 x 5 matrix with entries that are functions of the bgs V, while b is defined by 



(4.2.17). For instance, 



(4.2.20) 



A°(V) 



fu° 











\ 





u° 


QU 1 /!! 


QU 2 /U° 


QU 3 /U° 





n oi 


(R+P)U° 











n 02 





(R + P)U° 





\o 


n 03 








(R + P)U°) 



and similarly for the A fe (V), for k = 1,2,3. 

Remark 4.2.1. We reserve the use of matrix notation for the "upper half for two reasons. The first is 
that the "lower half involves constant coefficient differential operators, so when differentiating the "lower 
half equations, we don't have to worry about commutator terms, which are easily expressed using matrix 



notation as in (7.2.221, arising from differential operators acting on the coefficients. The second reason is 
that in future work, we plan to study the "lower-half in its original form as an inhomogeneous Klein-Gordon 
equation, but we will still use matrix notation for the "upper-half." 

Remark 4.2.2. A calculation gives that det (A°(V)) = -Q(R + P) 2 (U ) 3 , and in the Cauchy problem 
studied below, this formula will ensure that A° is invcrtiblc. 



5. The Geometry of the EN k System 

In this section, we discuss the geometry of the characteristics of the EN K system and relate the geometry 
to the speeds of propagation. 

5.1. The symbol and the characteristic subset of T*Ai. The symbol of the equations of variation 
at a given covector £ S T*Ai, the cotangent space of M. at x, is a linear operator on the space of variations 
V. This operator is obtained by making the replacements d\U — > on the left-hand side of the system 



(4.2.11 - (4.2.6). Here, U stands for any of the unknowns. The characteristic subset of the cotangent space 
at x is defined to be the set of all covectors £ <E T*A4 such that <j£ has a nontrivial null space. Thus, £ lies 
in the characteristic subset of T*A4 iff the following algebraic system has non-zero solutions V C K 10 : 



(5.1.1) 






= 






(5.1.2) 






= 






(5.1.3) 


(R + 




= 


Ci = 


1,2,3) 


(5.1.4) 






= 






(5.1.5) 






= 


(i = 


1,2,3) 


(5.1.6) 






= 0. 







The determinant of the linear operator o~£ at x, known as the characteristic form of the EOV and denoted 
by Q{x;£), is given by 



(5.1.7) 



Q(x;0 = (£o) 3 (^ a £a) VVT^^ 



where h 1 is the reciprocal acoustical metric, a non-degenerate quadratic form on T*Ai defined by 



(5.1.8) 
(5.1.9) 



(ft- 1 )"" = -a- 2 U»U u 
a = §(e'^P,Tf), 



1) [W, 



and the function § is defined by (3.3.81. The characteristic subset of T* is therefore equal to the level set 
(5.1.10) {UT:M\Q{x-i) = Q}. 
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Consequently, £ is an element of the characteristic subset of T*M, iff one of the following four conditions 
holds: 



(5.1.11) 
(5.1.12) 
(5.1.13) 
(5.1.14) 



W = 

(K i r v ^u = o 

£o = 0. 



Condition (5.1.111 defines a plane P* ~ in T*M, while conditions (5.1.121 and (5.1.131 define cones 



C* 



x,s{ound 



and C. 



x,U 



x,l(ight) ' 



respectively, in T*A4. Condition (5.1.141 also defines a plane P* in T*A4, and its 



presence is a consequence of our choice of dt4> as a state-space variable in our rewriting of the linear Klein- 



Gordon equation as a first order system. We refer to (5.1.11 1 - (5.1.141 as the four sheets of the characteristic 



subset of T*M. Fig. [T] illustrates the characteristic subset of T*M. In the illustration, we masquerade as if 
the domain of solutions to the EOV is R 1+2 , with the vertical direction representing positive values of £o- 




5.2. Characteristic surfaces and the characteristic subset of T X A4. A C 1 surface 5cM that is given 
as a level set of a function <!> is said to be a characteristic surface if at each point x S S, the covector £ with 
components £„ = <9„<I> for v = 0, 1, 2, 3, is an element of the characteristic subset of T*Ai. It is well-known 
(consult e.g. [9]) that jump discontinuities in weak solutions can occur across characteristic surfaces, and 
that characteristic surfaces play a role in determining a domain of influence of a region of spacctimc. 

There is an alternative characterization of characteristic surfaces in terms of the duals of the sheets P* - , 
Pg- 0) si an d C* v The notion of duality we refer to is as follows (consult e.g. [3]): To each covector £ 
in the characteristic subset of T*M. there corresponds the null space of £, which we denote by N%. This 
3-dimensional plane is a subset of T X A4, the tangent space of M. at x, and is described in coordinates 
as N% = {X £ TxMlZpX 11 = 0}. We define the dual to a sheet of the characteristic subset of T*M to 
be the envelope in T X M. generated by the JVj as £ varies over the sheet. The characteristic subset of the 
tangent space at x is defined to be the union of the duals to the sheets (5.1.111 - (5.1.141. A calculation 



of the envelopes implies that the respective duals to (5.1.111, (5.1.12), (5.1.131, and (5.1.141 are the sets of 



X 6 T X M such that in our fixed rectangular coordinate system (see Section 2.1 ) 



(5.2.1) 


X 


= XU for some AeR 


(5.2.2) 


K v x»x v 


= 


(5.2.3) 


9^X»X V 


= 


(5.2.4) 


X 


= A(l, 0,0,0) for some A e 


where 







(5.2.5) 



h^v = f fl> + (1 - o- 2 )U p U v 
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is the acoustical metric, a non-degenerate quadratic form on T X M.. The dual to P* ~, given by (5.2.1 1, is the 



linear span of U , and the dual to the plane P* , given by (5.2.4), is the linear span of ( 1 , , , ) . The dual to 
C* , given by (I5.2.2I) and labeled as C x s , is the sound cone in T X M, while the dual to C* ; , given by (15.2.31 



and labeled as C x j, is the light cone in T X A4. We refer to these subsets of T X M. as the four sheets of the 



characteristic subset of the T X M (noting that the degenerate cases (5.2.11 and (5.2.41 are lines rather than 
"sheets"). See Fig. [2] for the picture in M 1+2 , where the vertical direction represents positive values of X°. 



span (1,0,0,0) span U 




FIGURE 2. The Characteristic Subset of T X M 

It follows from the above description that for each £ belonging to a fixed sheet of the characteristic subset 
of T*Ai, N% is tangent to the corresponding sheet of the characteristic subset of T X A4. Therefore, we may 
equivalently define a characteristic surface as a C 1 surface S such that the tangent plane at each of its points 
x is tangent to any of the four sheets of the characteristic subset of T X A4. 



Remark 5.2.1. Note that C x , s lies inside C x j, but C* j lies inside C. 



x,s' 



5.3. Inner characteristic core, strict hyperbolicity, spacelike surfaces. The inner characteristic core 
of the cotangent space at x, denoted I*, is the subset of T*M. lying strictly inside the innermost sheet C* l . 
I* comprises two components, and we refer to the component such that each co-vector £ belonging to it has 
£0 > as the •positive component, denoted by 1* + : 

(5.3.1) T+ d = f U e T*M\^ < and £0 > 0}. 

A co-vector £ £ T*M is said to be hyperbolic for Q at x iff for any co-vector v not parallel to £, 



Q(x; A£ + v) — has real roots in A, where Q is given in (5.1.71. The set of hyperbolic co-vectors at x is 
equal to 2*; see Fig. [I] A co- vector £ S T*M is said to be strictly hyperbolid^for Q at x iff for any co- vector 
v not parallel to £, Q(a;; A£ + u) = has distinct real roots in A. As mentioned in Section [I] the EOV (and 
hence the EN K system) are (is) not strictly hyperbolic because of the repeated factors in the expression 



(5.1.71 for Q(x; ■), and because two of the sheets of the characteristic subset of T*M intersect. 

A C 1 surface S C M is said to be spacelike (with respect to the light cones C* ; ) if at each x & S, there is 
a co- vector £ belonging to 1* such that the tangent plane to S at x is equal to iVj. Based on the discussion 
above, it follows that S is spacelike at x iff the tangent plane to S at x is the null space of a co-vector £ that 
is hyperbolic for Q at x. 

5.4. Speeds of propagation. It is well-known that for first order symmetric hyperbolic systems, the speeds 
of propagation are locally governed by the characteristic subsets. For example, in the case that the charac- 
teristic subset of T*M at each x includes an innermost sheet, the domain of influence of a spacetime point 
x' is contained in the interior of the forward conoid in M. traced out by the set of all curves emanating 
from x' and remaining tangent to the sheets of the characteristic subsets of the T X M. that are dual to the 
innermost sheets of the characteristic subsets of the T*M. as the curve parameter varies; consult [21] for a 
detailed discussion of this fact. 

We will later illustrate the occurrence of similar phenomena in the EN K system. In this case, the innermost 
sheet at x is C* t , the dual of which is G x .i, the light cone in T X M.. Therefore, the forward conoid emanating 
from a spacetime point x' is the forward light cone in M. with vertex at x' . Thus, one would expect that the 



J For PDEs derivable from a Lagrangian, the notions of hyperbolicity, characteristic subsets, etc., have been generalized by 
Christodoulou [B] in a manner that allows one to handle characteristic forms that feature multiple roots. 
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fastest speed of propagation in the EN K system is the speed of light. This claim is given rigorous meaning 



below in the uniqueness argument (see Section 7.3.11 which shows, for example, that a solution that is 
constant in the Euclidean sphere of radius r centered at the point s £ M 3 at t = remains constant in the 
Euclidean sphere of radius r — t centered at s at time t > 0; see Remark |7.3.2| 

We contrast this to the case of the special-relativistic Euler equations without gravitational interaction, 
in which there is no Klein-Gordon equation governing the propagation of gravitational waves at the speed 
of light, and the set C* z does not belong to the characteristic subset of T*M. The inner sheet at x in this 
case is C* s , the dual of which is C x>s , the sound cone in T X A4, and the methods applied below can be used 
to show that the fastest local speed of propagation is dictated by the sound cones C x>s . This case is studied 
in detail in [SJ and [7]. 

5.5. Energy currents. The role of energy currents in the well-posedness proof is to replace the energy 
principle available for symmetric hyperbolic systems. After providing the definition of an energy current, 



we illustrate its two key properties, namely that it has the positivity property (5.5.21 below, and that its 
divergence is lower order in the variation V. 



5.5.1. The definit ion of an energy current. Given a variation V : M. 
defined in Section 



)10 



and a bgs V : A4 - 

we define the energy current to be the vectorfield J with components J , J- 7 , j 
in the global rectangular coordinate system given by 



aio 



4.2 



as 
2,3, 



JO clef ^2 



P- 



(4>i) 2 - 



P)U C 



u k u k 



(U k U 



k\2 , 



(^ 2 ) 2 + ftfcO 



(5.5.1) 



Q 



p2 



2WP+(R+P)W 



U k Uk 



(U k U k ) 2 



(t/°) 2 



Notation. In an effort to avoid cluttering the notation, we sometimes suppress the direct dependence of J 
on V and V and instead emphasize the indirect dependence of J on (t,s) through V and V by writing 

Terminology: We say that J is the energy current for the variation V with coefficients defined by the bgs V. 

Remark 5.5.1. The theory of hyperbolic PDEs derivable from a Lagrangian, and in particular the derivation 
of energy currents, is developed by Christodoulou in [BJ. For readers interested in studying Christodoulou's 
techniques, we remark that the Lagrangian density for (4.1.181 - ( 4.1.20[ ) (the first 5 scalar equations of 
the EN K system) is expressed in the original variables as pe 4 *. The energy current (5.5.11 is the sum of 
an energy current for the linear Klein-Gordon equation, which supplies the terms involving {(j)) 2 and (Vv) 2 , 
and an energy current used by Christodoulou in [7 to study the special-relativistic Euler equations without 
gravitational interaction. 

5.5.2. The positive definiteness of £ M J M for P > and £ € T* + . Given an energy current as defined by 



(5.5.1 1 and a co-vector £ e T*A4, the quantity £ M J M may be viewed as a quadratic form in the variations V 
with coefficients defined by the bgs V. We emphasize this quadratic dependence on the variations by writing 
£ M J M (V, V). One of the two key features of the energy current is that P > and £ 6 I* + together imply 
that the form ^J^V, V) is positive definite in V : 

(5.5.2) £ M >(V, V) > if £ e {C € T*(M) \ C^C < and Co > 0} and V + 0. 

A direct verification of this fact can be carried out, for example, by calculating the eigenvalues of the matrix 

of the quadratic form £ M J^(V,V). The eigenvalues depend on £ and are positive whenever P > and 



£ € T* + . As we shall soon see, inequality (5.5.21 will allow us to use the form £/ 1 .7'*(V, V) to estimate the 
L 2 norms of the variations, provided that we estimate the bgs V. 

Remark 5.5.2. Although later in this article we make use of the fact that V is a solution to the EOV, the 



inequality in (5.5.2) does not rely on this fact; it is an algebraic statement about ^^(V, V) viewed as a 
quadratic form on IR 10 . 
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5.5.3. The divergence of the energy current. If the variations V are solutions of the EOV, then we can 



compute duJ 1 ^ and use the equations (4.2.11 
derivatives of V : 



(4.2.6) for substitution to eliminate the terms containing the 



(5.5.3) 




2t)/ + 2 



Pg 
Q 



2t/ fc /i (fe) 



(C/°) 2 
U j h^U k U k 



jjk p 



2U k U k {R 




(IP) 



0)2 



That the right-hand side of (5.5.31 does not contain any derivatives of the variations is the second key 



property announced at the beginning of Section 5.5 



Remark 5.5.3. Given a spatial derivative multi- index a and an energy current J as defined in (5.5.1) such 
that the variation V is a solution of (4.2.1) - (4.2.10) with inhomogeneous terms (b, 1), where b and 1 are 



defined by (4.2.171 and (4.2.181 respectively, we define the higher-order energy current J& to be the energy 



current for the variation <9sV with coefficients defined by the same bgs V. The variations d&V are solutions 



of (4.2.1 1 - (4.2.101 with inhomogeneous terms (hs,dsl), where hs is defined in terms of b below through 



(7.2.21). Consequently, the expression for d^Ji is given by taking the formula (5.5.3) for c^J 7 * and making 
the replacements V — > 83V and (b, 1) — > (hs,dsl). 

6. Assumptions on the Initial Data 

We now describe a class of initial data to which the energy methods for showing well-posedness can be 
applied. The Cauchy surface we consider is {(^,s) S M. \ t = 0}. 



6.1 

V(s) 



An H 

def 



N 



perturbation of a quiet fluid. The initial data for the EN K system are denoted by V = 



(77, PjU 1 , • • • , ^3), where tpj — dj<j> for j 



1, 2, 3. We assume that the initial data V for the EN K sys- 
tem are constructed from initial data (77, j3, u 1 , ■ ■ ■ , ^3) in the original state-space variables (r],p, it 1 , • • • , ?/> 3 ) 
according to the substitutions ( 4.1.4[ ), (4.1.51, and (4.1.6). Additionally, we assume that outside of the unit 
ball centered at the origin in the Cauchy surface 



(6.1.1) V = V 

where <fi is the unique solution to 



(77, P, 0,0, 0,0, 0,0, 0,0), 



(6.1.2) 



K 2 + e 4 ^(3l(p,77)-3p) = 0, 



77 and p are positive constants denoting the initial entropy and pressure of the fluid outside of the unit ball, 
P = f e 4 ^p, and the function H is defined in (3.3.9). An initial state of this form is a perturbation of an 



infinitely extended quiet fluid, such that the perturbation is initially contained in the unit ball. Here we need 
the cosmological constant k 2 > in order to ensure that the EN K system has non-zero constant solutions of 
the form V. 

Because the standard energy methods require that the initial data belong to a Sobolev space of high 
enough order, we assume that 

(6.1.3) H77 - f]\\ H N + \\p-p\\ H N + \\u k \\ H N + \\(j)- 4>\\ h n+i + \\ipo\\ H N < OO, 

where N € N satisfies 



(6.1.4) 



N > 3. 



Note that (|6.1.3) implies that ||^j||ijN < 00 (j = 1,2,3). By Proposition 



from (6.1.3| that 
(6.1.5) 



A.6 



and Remark 



A.5 



it follows 



V|| ff jy < 00. 
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Remark 6.1.1. It is not necessary to assume that the initial deviation from the constant state has compact 
support. It is sufficient to consider initial data V that differ from V by a perturbation belonging to H N , 
such that that V(M 3 ) is contained in a compact subset of 0, where N is given by (6.1.41 and is defined 
in Section |6.2| We make the compactness assumption because it is useful for illustrating the speeds of 



propagation as discussed in Section 5.4 and because we plan to make use of this setup in future work. 

6.2. The admissible subset of state space and the uniform positive deflniteness of J°. In this 
section we discuss a further positivity restriction that we place on the initial data. We will see in Section 



7.2.3 that this positivity condition is propagated for short times during an iterative construction of solutions 
to the linearized EN K system. Since it plays a key role in our future analysis, we discuss here the implications 
of this positivity restriction regarding the uniform positive deflniteness of the energy current, viewed as a 
quadratic form in the variations. 

6.2.1. The definition of the admissible subset of state-space. In order to avoid studying the free boundary 
problem and in order to avoid singularities in the energy current, we assume that the initial pressure, energy 
density, and speed of sound are uniformly bounded from below by a positive constant. According to our 



assumptions (3.3.51 on the equation of state, to satisfy these requirements, it is sufficient to consider initial 
data for the EN K system such that V(IR 3 ) is contained in a compact subset of the following open subset 
of the state-space R 10 , the admissible subset of state-space: 



(6.2.1) 



= {V £ R w \rj > Q,P > 0}. 



We therefore assume that V £ 0\ and V(R 3 ) C Ox, where Ox is a precompact open set with Ox <s O. We 
then fix a precompact open subset 02 with conve^ closure satisfying 0\ <s O2 C 02 <e 0\ our goal is to 
show the existence of a solution that remains in 2 for short times. 

6.2.2. The uniform positive deflniteness of J° . Most of the technical exposition below is devoted to obtain- 
ing control over ||V(f)||jjJv, where V is a solution to the EOV defined by a bgs V. Instead of trying to 
estimate |jV(t)|| £ 2 directly, it is advantageous to estimate || J°(i)||i,i, where J is an energy current for V 
with coefficients defined by the bgs V, since the divergence of J is lower order in V. A similar remark applies 
to estimating ||<95V|| L 2 using higher-order energy currents J$. We shall see that || J°(^)||li can be used to 
estimate ||V(t)||? 2 from above and below provided that J° is uniformly positive definite independent of the 
bgs V. More precisely, we claim that there exists a Cq 2 with < Cq 2 < 1 such that for any variation V 
and any bgs V contained in 02, we have 

|V| 2 < j°(V,V) < 1 



(6.2.2) 



Q) 2 



c©2 



IVI 



To prove ( |6.2.2[ ), recall that J is defined by fl5.5.1| and note that (1,0,0,0) £ T*+ by ( |5.5.2| . The uniform 
continuity of J (which we momentarily view as a function of (V, V)) on the compact set 02 x {|V| = 1} 



implies that there exists a Cq 2 with < C<g 2 < 1 such that (6.2.2 1 holds whenever V(t, s) £ O2 and |V| = 1 



Since the inequalities in (6.2.21 are invariant under any rescaling of V, it follows that we may remove the 
restriction IVI 



1. 



7. The Well-Posedness Theorems 



In this section, we state and indicate how to prove our two main theorems. We have separated the proof 
of well-posedness into two theorems since the techniques used in proving each are different. Statements of 
the technical estimates involving the Sobolev-Moser calculus have been placed in the Appendix so as to not 
interrupt the flow of the main argument. 



Proposition A.8 requires the convexity of 02- Without loss of generality, we may choose it to be a cube. 
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Theorem 1. (Local Existence and Uniqueness) Let V(s) be initial data for the EN K system (4.1.181 - 
(4.1.26) that are subject to the conditions described in Section^ Then there exists a T > such that (4.1.181 
- (4.1.261 has a unique classical solution V(t, s) on [0,T] x R 3 satisfying V(0,s) = V(s). The solution is of 
the form V = (rj, P, U 1 , U 2 , U 3 , 4>, do<j>, d\(f>, 82(f), d^cj)) and satisfies 
V([0,T] x R 3 ) C 2 . Furthermore, 

V S Cl{[0,T] x R 3 ) n C°([Q,T],H$) n C^^T],^" 1 ), and consequently 

e c 2 h ([o,T] x M 3 ) n C°([0,T],iff +1 ) n tfHM,^ ) n c^M,-^ -1 )- 

Proof. As discussed in Section [7~T] our abbreviated proof of Theorem [l] is located in Section 7.2 



□ 



Remark 7.0.1. In the discussion below, we sometimes denote the solution from Theorem [T] by ~V S ol for 
clarity. 

Corollary 7.0.1. The interval of existence [0,T] supplied by the Theorem^ depends only on the set O2 
from Section ||(°)V||jy.iv+i, and the constant A chosen in (7.2.71 - ( 7.2.9[ ) below. Here, (°)V denotes the 

mollified initial data as described in Section 



7.2 



Furthermore, the set O2, the mollified initial data C'V, 
and constant A can be chosen to be independent of all initial data varying in a small H N neighborhood o/V. 

Therefore, if we define _B a (V) = {V S | ||V — V||jjn < y}, then there exist S > andT' > (depending 

onW) such that any initial data V belonging to Bg(V) launch a unique classical solution V that exists on 
the common time interval [0, T'] and that has the property V([0,T'] x M 3 ) C 02- 



Proof. The corollary follows from the proof of Theorem]!] See in particular Remark 7.2.1 and Remark 7.2.2 
below. □ 

Corollary 7.0.2. The norms \\\ V |||^fiy T and \\\ d{V \\\h n - 1 ,t of the solution from Theorem^depend only 
2 , || (0) V|| H iv + i, and A. Furthe rmore, there exists a K > such that any initial data V belonging to the set 
B$(V) defined in Corollary 7.0.1 launch a unique solution V that satisfies the uniform bound 

d t V ||| ff «-i^</Wa 2 ,|| (0) V|| H « +1 ,A, ( 5), 



(7.0.3) 



V 



where T and T' are as in Corollary 7. 0. 1 



7.2.1 



and 



Proof. The estimates for ||| V \\\ h n t and ||| V \\\ h n t , follow from Corollary 7.0.1 Proposition 

the fact that the sequence of iterates {'• m ^V(i)} constructed below converges strongly in L 2 ^ and weakly 
in H^. to V(i); consult [55] for the missing details. We then use the EN K equations to solve for the time 
derivatives together with Proposition A. 6 and Remark A. 5 to obtain the estimates for 
II 9tV \\\h n '- 1 ,t an d HI <9(V Hljyjv-i x> ■ D 

Theorem 2. (Continuous Dependence on Initial Data) Let V(s) be initial data for the EN K system 
(4.1.18) - (4.1.261 that are subject to the conditions described in Section^ and l et V b e the solution existing 
on the time interval [0,T] furnished by Theorem^ Let Bs(V) be as in Corollary 7.0.1 Let {V™ 1 } C Bs be a 
sequence of initial data withlimm^aa ||V m — V||^fw = 0, and let~V m denote the solution to ( |4.1.18 1 - (4.1.261 
launched by V m . Then for all large m, the solutions V™ 1 exist on [0, T], and linim^oo ||| V m — V |||jjn t= 0. 



Proof. Our proof of Theorem [2] is located in Section 7_A 



a 



Remark 7.0.2. It is unknown to the author whether or not the continuity statement from Theorem [2] 
can be strengthened to one of Lipschitz continuity or Holder continuity. However, using Burger's equation 
dtu+ud x u = 0, Kato [TB] provides a counterexample in which the map from the initial data uo G H a (M.) to the 
solution u(t) £ C([0,T], H a ) is not Holder continuous with any positive exponent; such a counterexample is 



explicitly constructed for a > 2. On the other hand, inequality (7.3.271 below shows that for the EN K system, 
the map from the initial data to the solution is a Lipschitz-continuous map from into C([0, T], H^." 1 ). 
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7.1. A discussion of the structure of the proof of the theorems. We prove local existence by following 
a standard method described in detail in Majda's book [22]: we construct a sequence of iterates {( m )V(t, s)} 
that converges to the solution V so /(i, s). To construct the iterates, we first define a sequence of C°° initial data 
{WV} such that ( m )V(M 3 ) <s 2 and lim ro ^oo II (m) V - V \\ H x= 0. The advantage of smoothing the data 
is that all of the iterates are C°° , thus allowing us to work with classical derivatives during the approximation 
process. Then beginning with (°)V(i,s) = f (°)V(s), we inductively define ( m+1 )V(i,s) as the unique solution 
to the linearization of the EN K system around ( m )V(t,s) with initial data ( m+1 )V(0,s) = ( m+1 )V(s). As a 
consequence of the theory of lineaiQ PDEs (consult [S]), each iterate (™ l )V is known to possess a classical 
solution on a strip [0,T m ] x M 3 , on which it satisfies, for every real number N', £ C°([0, T m ], H$ ). 

Here, T m is any real number such that ( m-1 )V([0, T m ] x K 3 ) C C>2, which ensures that the sequence of proper 
energy densities is bounded from below by a uniform constant and therefore precludes singularities in energy 
the currents we use during the linearization process. 

In order for the limiting function "V so i to be defined on a strip, it is obviously necessary that we show 
that the sequence of time values {T m } can be bounded from below by a positive constant J"*. To this end, 
we examine the EOV satisfied by ( m ) V — (°) V and its partial derivatives, and we control the growth in 
of HI ( m )V — (°)V |||.h-jv,t. uniformly in to using energy currents. According to the above paragraph and 
the Sobolev embedding result H 2 (R 3 ) C C£(M 3 ), it follows that if ||| WV - \\\ H ",t, is small enough, 
uniformly in to, then T* may be selected as a uniform lower boun d on the T m . Our detailed proof of the 
control of the terms ||| f m 'V - f°'V ||| Tt is given in Proposition 



7.2.1 



below and uses the Sobolev-Moser 

3 

2 - 



calculus inequalities, which are refined versions of the fact that for N' > §, H N (M 3 ) is a Banach algebra. 
Their purpose is to control the L 2 norms of terms of a product form, based on known Sobolev regularity 
of each factor in the product. We state the relevant Sobolev-Moser estimates in the Appendix and give 
references for readers interested in the proofs. 



Our proof of Proposition 7.2.1 illustrates the relevant techniques for obtaining Sobolev estimates from the 
method of energy currents. Instead of completing the existence proof, which requires arguments similar to 
the ones used in proving this proposition, we refer the reader to Majda's local existence proof for symmetric 
hyperbolic systems [55]; the only necessary modification to Majda's proof is to use the method of energy 
currents in p lace o f the energy principle for symmetric hyperbolic systems. 



In Section 7.3.1 we show uniqueness and H — Lipschitz-continuous dependence on the initial data. The 



methods used in this argument are similar to the methods used to prove Proposition 7.2.1 so we provide fewer 
details. We consider the EOV satisfied by the difference of two solutions V and V to the EN K system, and 
then use an appropriately defined energy current to bound the growth of ||| V — V |||jjw-i T by a constant 
times exponential growth in T. We show that the constant depends on the initial data and is bounded 
from above by another constant times ||V(0) — V(0)||#-jv-i , thus implying uniqueness and H 1 ^^ 1 — Lipschitz- 
continuous dependence on the initial data. Our abbreviated proof of Theorem [T] is complete at the end of 
this section. 

Our proof of Theorem [2] requires some machinery from the theory of evolution equations in a Banach 
space. The basic method is due to Kato [13], and most of the technical results we use in this section are 
merely quoted from his papers. We find it worthwhile to prove Theorem [2] because aside from Kato's work, 
we have had difficulty locating this result in the literature. 

7.2. An abbreviated proof of Theorem [lj As described in Section [7TT] we produce a sequence of iterates 
{( m )V(£,s)} that converges to the solution V so i(t,s). 

7.2.1. Smoothing the initial data. We begin by smoothing the initial data V, which we assume are of the 
form described in Section[6] so that we can work with classical derivatives. Let vE'(s) be a Friedrich's mollifier; 
i.e. * € C C °°(]R 3 ), supp(*) C {s| |s| < 1},* > 0, and / * d 3 s = 1. For e > 0, we set * e (s) = e _3 *(f) and 
define * e V e C°°(]R 3 ) by 

(7.2.1) * e V(s) d = / * e (s-s')V(s')dV. 



^The exposition on linear theory in [9] makes use of the symmetric hyperbolic setup to obtain energy estimates for the linear 
systems. We may obtain similar energy estimates for the linearized EN K equations by using energy currents of the form ( |5.5.1| ; 
the proof of Proposition |7.2.T] below illustrates the relevant techniques. 



19 



The following properties of such a mollification are well-known: 



(7.2.2) 
(7.2.3) 



lim m e y-y\\ H N 

3{e > A C(V) > 0} 3 < e < e 







|* £ V-V|| L2 <eC(V)||V|| ff i 



We will choose below an eo that is at least as small as the one in (7.2.31. Once chosen, for a given m G N, 
we define 



(7.2.4) 
(7.2.5) 
(7.2.6) 



e, n — I eo 

My = f * e V 

(™)W = * f W, 



where W denotes the first 5 components of V. 

By Sobolev embedding, by the assumptions on the initial data V, and by the mollification properties 



above, 3{A > A eo > 0} (at least as small as the eo in ( 7.2.3 1) such that 



(7.2.7) 
(7.2.8) 
(7.2.9) 



||V-<°>V| 
||V-<°>V| 

|(m)y _ (0)y| 



<r A 

< A ^ V(M 3 ) c & 2 
A 



H 



H' 



holds for 771 > 0, 



where C<g is defined in (6.2.21 



Remark 7.2.1. It is a standard result that if e > and N' is any real number, then ^ e V £ H^' (E 3 ). We will 
make use of this remark below, for in the local existence proof, we will need to differentiate the equations 



(7.2.15) - (7.2.201 N times and utilize Sobolev estimates; since several terms from these undifferentiated 



equations already contain one derivative of the smoothed initial data, our estimates will involve |p "^"llij-^ 



See e.g. (7.3.61 and (7.3.9) 



Remark 7.2.2. If we are considering initial data V in a small enough H N neighborhoo d JV o f the initial 
data V, we can use a fixed smoothed function (°)V in place of each (°)V in Proposition 7.2.1 below, and 



choose A to be uniform over the neighborhood. For what then enters into the proof of local existence for 
the initial data V are the quantities ||^°- ) V|| ff ]v+i and ||( m )V — ^V||#jv, and the latter quantity is easily 
controlled by the inequality 

l|M V - Wy\\ HN < ||(" l )v - y\\ HN + ||v - y\\ H N + \\v - (0) v|| H « ; 



(7.2.10) 



once we fix an appropriately chosen smoothed function (°)V and a corresponding A satisfying (7.2.71 and 



(7.2.8 1, we may independently adjust the mollification of each V belonging to Af so that the right-hand side 



of (7.2.101 is < Cq 2 A/2 for m > 0. This estimate would then enter into our proof in inequality (7.2.31 1. We 
also note that this remark is relevant for Corollary |7.0.1| above. 

7.2.2. Defining the iterates. Consider the iteration scheme described in Section [7.1| The components of the 
iterates are denoted by ( m )V = (y^T}, ( m )p, • • • , ^ m 'ips\ , and we use the notation ( m )W to denote the first 
five components of ( m 'V. Linear existence theory implies that each iterate ( m+1 )V is a well-defined, smooth 
function with ||( m+1 )V(i) - (0) V|j ff jv < oo for < t < T„ 
which HI ( m )V - ( °)V | \ \h n t < A holds. 



Here, by (7.2.81, [0, T m ] is any time interval for 



7.2.3. The uniform time estimate. As discussed in Section 7.1 we show the existence of a fixed > such 
that HI ( m )V — (°)V \ \\jjn Tt < A for 777 e N, thus ensuring that each iterate is defined for a uniform amount 
of time and remains inside of O^. We state a slightly stronger version of this result as a proposition: 
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Proposition 7.2.1. Let A denote the constant defined in (7.2.7 1 
L > such that each of the iterates ^ m ^V(i,s) satisfies 

(7.2.11a) I]] Mv-Wv 

(72.11b) |||5 t ( (m) V) 



(7.2.91. Then there exist T* > and 



Iff- 



satisfies (7.2.11a 



satisfies (7.2. 11a I and 



Proof. We proceed in our proof of Proposition 7.2.1 by induction on m, noting that (°)V(t,s) = (°>V(s) 
and ( 7.2. lib I with any T* > and any positive number L. We thus assume that ( m )V 
7.2.11b I without first specifying the values of or L. At the end of the proof, we 
will show that we can choose such a T* and an L, both independent of m, such that energy estimates imply 
the inductive step. To obtain the estimates stated in the proposition, it is convenient to work not with 
the iterates themselves, but with the difference between the iterate and the smoothed initial value. Thus, 

for each mSNwe define 



referring to the notation defined in (7.2.51 and (7.2.61 

(7.2.12) V(i,s) = (m+1) V(t,s) - 

(7.2.13) W(t, s) = ( m+1 >W(<, s) 

(7.2.14) V = (m) V. 

We have used the notation V and V suggestively: it follows from the the definition of the iterates, 



<°>V(s) 
- !°>W(8) 



definition ( |7.2.12[ ), and definition (|7_2_14j that V is a solution to the EOV ( |4.2.1[ ) - ( 4.2.10| | defined by 
the bgs V with initial data V(0,s) = ( m+1 'V(s) - (°)V(s). Our notation (|7.2.12| - (|7.2.14| is therefore 



consistent with our notation for the EOV introduced in Section 4.2 Recalling also the notation (4.2.171 



and (14.2.18 1 introduced in Section 4.2 the inhomogeneous terms in the EOV satisfied by V are given by 



(b,l) = (/,<?, 


• • , l^), where for j = 1, 2, 3, 




(7.2.15) 


f=-U k dk[Wfi\ 




(7.2.16) 


g=-U k d k [^P] - Q8 k [^U k ] + (4P- 




(7.2.17) 


h U) = -(R + p)u k d k [ w w] - n fc ^[ (0) p] 


+ (3P - R)0^% 


(7.2.18) 


l^ = n 2 ^+R-3P-d k [^ k ] 




(7.2.19) 


z« = dJ[(0)^ ] 




(7.2.20) 


= 





As explained in Section |4.2[ for each spatial derivative multi-index a with 
< \a\ < N, we may differentiate the EOV with inhomogeneous terms (b, 1) to which V is a solution, 
obtaining that 83V is also a solution to the EOV defined by the same bgs V with inhomogeneous terms 
(hs,dsl). The inhomogeneous terms b,; are given by 

(7.2.21) b s d = f A°d s {(A )- 1 ^ + k s , 

where 



(7.2.22) 



k^ = A 



(A°)- 1 A k d k (d s W) - d s ( y (A r 1 A k d k w) 



for < I a I < N. Note that we ha ve suppressed the dependence of the A v {-) on V 
As discussed in Section 



6.2.2 



we will use energy currents to control ||| V | ] | jyjv,x • We state here as a 
lemma an important differential inequality that allows us to proceed with our desired Sobolev estimates. Its 
proof is based on the key properties of energy currents described in Section [5.5| and the divergence theorem. 

Lemma 7.2.2. (See Fig. [Jj Suppose r > T > 0. For < t < T, let 

^t,r-t — { x <= A^|x° = t, x k x k < r — t} denote the Euclidean sphere of radius r — t centered at (t, 0, 0, 0) in 
the flat hyper surf ace {a; = i}, and let 

M t<r = {x G M\Q < x" < t, x k x k = r — x } denote the mantle of the past directed, truncated light cone with 



lower base S r and upper base T, tr _ t . Let V be a solution to the EOV (4.2.11 - (4.2.101 defined by the bgs 
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V, and assume that 

V([0, T] x R 3 ) C C>2- Let J be the energy current (5.5.1 ) for the variation V defined by the bgs V, and define 



£(t;r)= J / J°(t,s)d 3 s 

v s t , r _ t 

(7.2.23) 



1/2 



The 



2£(t;r)j t £(t;r)< j j"(t, s) d 3 s. 



Remark 7.2.3. We note that our use of i in the statement of Lemma [7.2.2| as a constant value taken on by 
the generic spacetime coordinate x a is inconsistent with our usual notational convention for spacetime points 
defined in Section |2.1[ in which t and xq are both used in the same manner as generic coordinate variables. 



Proof. By the divergence theorem, we have that 



£ 2 (t;r)= J j°(t,s)d 3 s = jj°(0 : s)d 3 s 



(7.2.24) 



(n(x),J{x)) E dH{x) + I II d^j f *(t',s)d 3 s)dt'. 



Aft,. 



t'=t 
t'=0 




Figure 3. The Surfaces of Integration in Lemma [7. 2. 2| 



Here, n(x) is the Euclidean outer normal at x £ M t . r to the mantle of truncated cone, (h(x), J(x)) e 
denotes the Euclidean inner product of n(x) and J(x) as vectors in R 4 , and Tt is the Hausdorff measure on 
the mantle of the cone. For each normal vector n(x), let denote the co-vector belonging to T*A4 such 
that n(x)£(X) = fitx)^^^ = ( n ( x ),X)E holds for every X E T X A4. By the positivity condition (5.5.2), co- 
vectors £ belonging to T* + satisfy ^ M J M (V, V) > for all non-zero variations V. Since for each x £ M t r , the 
co- vector fi(x)£ belongs to the boundary of T* + , which is the positive component of the cone C* ; , continuity in 
the variable £ implies that (n(x), J(x))e — n(x)^J tl (^ |V) > holds for x £ M t r . Furthermore, if t\ < 



then M tl . r C M t2 . r . From these facts it follows that 



- J (n(x), J{x))e d7i(x) is a decreasing function of t 

M t , r 



on [0, T]. Lemma 7.2.2 now follows from differentiating each side of (7.2.241 with respect to t and accounting 
for this decreasing term. Fig. [3] illustrates the setup in K 1+2 , where the vertical direction represents positive 
values of t. □ 



Returning to the proof of the proposition and recalling that we are using definitions (7.2.12 1 and (7.2.141 



to define V and V, we let J3 denote the energy current for the variation 9,;V defined by the bgs V. For 
notational convenience, we allow a to take on the value 0, in which case Jg is defined to be the energy 
current in the variation V defined by the bgs V. 
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As in Lemma 7.2.2 we define for any T» > and r > T* the following functions of t on [0, T*] : 

def / r ... . -o \ 1/2 



(7.2.25) 



£ S (i;r)=(^ 1 j° s (t,s)d 3 s 
S(t;r;iV) d = f ( £ £|(t; 



0<|<3|<7V 



E / JSM^ 

o<|5|<w E ; 



Then with Cg 2 defined in (6.2.21, we have that 
(7.2.26) 



C 02 E'(t;r;N) < ||V(t)||^ (E , < C£ (t; r; N) . 



Additionally, by Lemma 7.2.2 we have the following inequality for < t < 



(7.2.27) 



2E(t;r;N)-E(t;r;N) < E 



0<\3\<N^' t 



cUje(t,s) d's 



The technically cumbersome aspect of the proof of Proposition |7.2.1| is bounding the right-hand side of 
(7.2.27) by a constant times E(t;r;N) + E 2 (t;r;N), which then allows us to use Gronwall's inequality to 
exponentially bound from above the growth of E(t; r; N) in t. We prove some of the technical points in 
lemmas 7.3. 1| and 7.3.2 below, so as to not disrupt the main argument. The keys to proofs of lemmas [7.3.1 
and 7.3.2 are Sobolev-Moser calculus inequalities, special versions of which are stated in the Appendix. In 
the following argument, C = C(N,6 2 , || (0) V|| h n+i , A, L), even though we have not yet chosen L. By Lemma 
7.3.1| we have that 



E / d M (j£{t,B)) d 3 s<C- [||V(i)|| HN(Str _ t) + ||V(t)|| 2 ffN(Sfi ._ t) ] 



o<|S|<at e ; 



(7.2.28) 



< C ■ [C 6 l /2 E(t; r; N) + C^E 2 (t; r; N)} , 



where in the second inequality we have used ( 7.2.26[ ). Combining (7.2.271 with (7.2.281, and applying 
Gronwall's inequality, we have for < t < that 



(7.2.29) 



E(t-r-N) < [E(0;r;N) + C-(2C 62 )- 1/2 t] ■ [exp(C • (2C 6a )-H)] , 



and consequently by (7.2.26), that 

(7-2.30) l|V(t)|| ff » (St _ t) < C l[\\V(O)\\ HN ^ Or) +Ct] ■ [exp(Ci)]. 

Letting r — ► oo, taking the sup over t € [0,T»], and using ( |7.2.9 1, we have that 



III V \\\ h »,t.< C0l[\\V(O)\\ H » + CT m ] •expfCT.) 
(7.2.31) < [A/2 + C(iV,O2,|| (0) V|| H ^ + i,A,i)-T t ] ■exp(C(N,6 2 ,\\ {0) V\\ H N + i,A,L)-T*) 

To make a viable choice of L, 
(7.2.32) 



we first assume that right-hand side of (7.2.31 1 is < A 



which implies the inductive step (7.2.11a) for ( m+1 )V. Using assumption (7.2.321 as a hypothesis, Lemma 



7.3.2 



implies that there exists an L{N,0 2 , ||^V||#n, A) > such that 
(7-2.33) HI <9 t ( (m+1) V) |||jy*-i,T„< L. 

For this fixed choice of L, we can implicitly solve for aT, > such that the right-hand side of inequality 



the induction argument. This completes the proof of Proposition |7. 2.1] 



(|7.2.31|) is in fact < A, thus justifying the assumption (|7.2.32|) and the conclusion (7.2.331, thereby closing 

□ 
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7.3. Proofs of Lemma 7.3.1 and Lemma 7.3.2 We now state and prove the two technical lemmas 



quoted in the proof of Proposition |7.2.1 



Lemma 7.3.1. Assume the hypotheses and notation of Proposition 7.2.1 In addition, assume that 
\\d t V(t)\\ H N-i < L. Then 



(7.3.1) 



d^(t)\\ LH ^ r _ t) < C(N,d 2 , ||(°)V|| <+I ,A,L) • (||V(f)|| H * + ||V(t)|| 2 ^ 



Proof. We use here the definitions (7.2.12 I and (7.2.14) from Proposition 



7.2.1 



Recall that <9 3 V is a solution 



to the EOV defined by the bgs V with inhomogeneous terms (b^c^l), and that J$ is the energy current 
for ds V defined by the bgs V. Furthermore, 



(7.3.2) 



!|V- (0) V|| H iv < A 



holds by the induction ass umptio n from the proposition 
By (5.5.31 and Remark 



5.5.3 



the expression for d^J'i consists of terms that are either precisely linear 
or precisely quadratic in the components of the variation <9sV. The coefficients of the quadratic variation 
terms are smooth functions with arguments V and D"V . Examining the particular form of these coefficients 
and using the fact that V([0,T] x R 3 ) C 6 2 , we see that their L°° norm is bounded by C{6 2 )\\DV\\ L °°. 
By assumption, ||V - (0) V|| H iv < A and \\d t V{t)\\ H N-x < L. Therefore, by Sobolev embedding, ||DV||l« < 
C(N,\\WV\\ h n,A,L). These facts imply that the L l (E iir - t ) 

norm of the terms involving the quadratic 

variations is bounded from above by C(N,0 2 , W^VWh? , A, L)\\dsV\\ 2 L2 ^ t \. 

The coefficients of the linear variation terms are linear combinations of products of the components of 
(b<s, dsl), where h& is defined in (7.2.21 ), with smooth functions, the arguments of which are the components 
of V. Since 

V([0,T] x R 3 ) c 6 2 , the smooth functions of V arc bounded in L°° by C{6 2 ). Therefore, by the Cauchy- 
Schwarz integral inequality for L 2 , the L 1 (_Ht,r-t) norm of the terms depending linearly on the variations is 
bounded from above by C(£>2)||(b5,<95l)|| L 2||<9,3V|| L 2( Sf r t y To complete the proof of (7.3.11, it remains to 
show that for < \a\ < N, we have that 



(7.3.3) 



(b 3 ,5 5 l)|| L2 < C(N,<D 2 , \\Wv\\ B „ + i,A) ■ (1 + ||W||^ 



The proof of (7.3.31 will follow easily from the propositions given in the Appendix. 

Concerning ourselves with the ||b<j||/,2 estimate first, we claim that the term A°d s ((A ) _1 b) from ( |7.2.21| 
satisfies 



(7.3.4) 



\A°d s ((A )- 1 b)\\ L2 <C(N,6 2: \\^V\\ H , +1 ,A). 



We repeat for clarity that b = (/, g, h 1 , h 2 , h 3 ), where the scalar- valued quantities /, g, h 1 , h 2 , h 3 are defined 
in (|7.2.15| - ( |7.2.17| . Since \\A°(V)\\ L ^ < C(Q 2 ), to prove ( |7.3.4[ ), it suffices to control the L 2 norm of 



ds ((A u ) b) . Using Proposition 



A.6 



and Remark 
that 



A.5 



with (A^) '"playing the role of F in the proposition 



and b playing the role of G, we have 
(7.3.5) WdsdA )- 1 ^ < IKAVbH^ < C(AT,0 2 , ||(0)v|| H ^,A)||b|U, 



Furthermore, Proposition |A.6| and Remark |A.5| imply that 

(7.3.6) \\b\\ HN <C(N,6 2 ,\\Wv\\ h n +1 ,A). 



Combining (7.3.5) with (7.3.6) proves (7.3.41. 
We next claim that the kg from (7.2.221 satisfy 

||k a || L 2 <C(iV,a 2 ,||( )V|| H «,A)||W|| 



(7.3.7) 
Ag 



Again, since ||A°(V)|| L °o < C(0 2 ), to prove ( |7.3.7| , it suffices to control the L 2 norm of (A°)~ 1 A k d k (d s W) - 
1 A k dk^V) ■ By Proposition 



A.9 



and Remark 



A.7 



with (A ) 1 A h playing the role of F in the 
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proposition, and <9fcW playing the role of G, we have that 
(7.3.8) \\(A°)~ 1 A k d k (d s W) - d s ((A°)~ 1 A k d k W)\\ L ^ 

<C(N,0 2 ,\\^V\\ h »,A)\\V^W\\h»-i, 



from which (7.3.71 immediately follows. 

To finish the proof of (7.3.31, we will show that 

(7.3.9) ll^llff" <C(N,d 2 ,\\^V\\ H N +l ,A) (z = 0,1,2,3,4). 

For 1^,1^1^,1^ defined in ( |7.2.19| ) and ( |7.2.20| ), the claim is trivial. To estimate the component ; (0) , 
defined in (7.2.181, we first rewrite 

(7.3.10) /(°) = K 2 (0-, 



(R-R)-3(P-P)~d k [(°H k ], 



where P = e^p and R = e 4 ^Ji(e 4 ^P,r]), the function 3? is defined in (4.1.51, and p and r] are constants 
defined in Section [Tj] In equation (7.3.10), we have made use of (6.1.21, which is the assumption that 
k 2 4> + R-'SP = 0. Since 

(7.3.11) +3\\(P - P)\\h» + \\d k { m ^k}\\ H * <C(||(°>V|| <+1 ,A), 
we only need to show that 

(7.3.12) \\R-R\\h» <C(N,6 3 ,\\^V\\ hS ,A). 



This follows immediately from definition ( 4.2.10[ ), Proposition A. 8 and Remark A. 6 



Inequality (7.3.3) now follows from combining (7.2.21), (7.3.4), (7.3.71, and (7.3.91; this completes the 



proof of ( 7.3.1 1. 



□ 



Lemma 7.3.2. Assume the hypotheses and notati on of Proposition 7.2.1 Also assume the induction hypoth 



esis HI ( m )V — (°)V |||i/»,T»^ A from Proposition 
Then 



7.2.1 



Assume further that ||| ( m+1 IV-(°'V \\\ h n T < A. 



(7.3.13) 



dt ( (m+1 *>V) | \\ H s- liTm < L(N, 2 , || (°) V\\ HN , A) . 



Proof. By Remark 4.2.2 we may solve for 9t(' m+1 )W) : 

(7.3.14) <9 t ( (m+1) W) = (A°(( m )V)) _1 (b - A fe (( m >V)0 fc (( m+1 >W)), 

where the function b denotes the inhomogeneous terms from the linearized EN K equations satisfied by 
(m+i) W; Le>) b=«8(( m )V), where 

(7.3.15) 58(-) = (5(-),«(')^ (1) (-)^ (2) (-)^ (3) (-)) 

is an array-valued function, the scalar-valued functions (£>, ■ • • , Sj^ are defined in (4.2.11 1 - (4.2.13 1, and 
the A»(-) are defined in ( |4.2.19 l. 

Using the hypotheses of the lemma, we apply Proposition |A.6| and Remark |A.5| to the right-hand side of 
(7.3.14), concluding that 



(7.3.16) 



9t ((-+i) W ) |||^-i )r .<L(JV ) <5 2 ,||P ) )V|| ffy ,A). 



Likewise, an argument similar to the one used to prove (7.3.91 gives that 

(7.3.17) |||9 t ((™ +1 ^( m+1 Vo,--- , (m+1) V> 3 ) \\\ h n- 1 ^<L{N,6 2 ,\\^V\\ hS ,A). 



Combining (7.3.161 and (7.3.171 proves (7.3.131. 



□ 
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7.3.1. Uniqueness and H 



N-l 



- Lipschitz-continuous dependence on initial data. We now prove the prove the 



uniqueness of solutions to the EN K system and show that the solution is an H N 1 — Lipschitz-continuous 
function of the initial data. Let V denote initial data that launch a solution V of the EN K system as 
furnished by the existence aspect of Theorem |T| Let S, Bg(V),T' , and K(N,0 2 , \\^'V\\ h n+i ) A,S) be as in 



corollaries 



7.0.1 



and 



7.0.2 



Assume that the initial data V belong to Bg , and let V be a s olution of the EN K 

We now define 



system with initial data V existing on the interval [0,T'] as furnished by Corollary 
(7.3.18) V= f V-V. 



7.0.1 



It follows from definition ( |7.3.18| that V is a solution to the EOV ( |4.2.1| - ( |4.2.10| defined by the bgs V 
with inhomogeneous terms given by (for j = 1, 2, 3) 



(7.3.19) 
(7.3.20) 

(7.3.21) 
(7.3.22) 
(7.3.23) 
(7.3.24) 
(7.3.25) 



f 
9 



h (j 



: (17" - U^O^ 

: (UK - U»)d M P+ [QU k /U° - QU k /U a ]d Q U k 
+ {Q- Q)d k U k + (4P - SQ)U»^ - (4P - 3Q)U^ 

= [(r + p)u" -(R + p^^df.w + (n™ - n w )a M p 

+ (3P - ^fi™^ - (3P - P)IP J > M 
;(°) = n 2 (4> - 4>) + (R - 3P) - (R - 3P) 
lV> = 



and we denote them using the abbreviated notation b and 1 defined in (4.2.171 and (4.2.181 



By combining Proposition A. 6 Remark A. 5 Proposition A. 8 and Remark A. 6 (noting the particular 



manner in which the inhomogeneous terms depend on the difference of functions of V and V) , we have that 
(7.3.26) HI (b,l) ||U*-i T '< C(N,(D 2 ,K) ||| V 



Without providing details, we reason as in our proof of Proposition 7.2.1 using (7.3.261 in place of (7.3.61 



and (7.3.9) to arrive at the following bound: 



(7.3.27) 



V 



I iff 



<C£||V(0)|| ffJ «r-x ■exp(C(N,d 2 ,K)-r) 



where V(0) = V - V. 

We now observe that ( 7.3.27 1 implies both the uniqueness statement in Theorem[l]and the P w_1 — Lipschitz- 
continuous dependence on the initial data mentioned in Remark 7.0.2| 



Remark 7.3.1. We cannot obtain an estimate analogous to (7.3.271 by using the norm in place of the 
i7^ _1 norm; the inhomogeneous terms (7.3.191 - (7.3.251 already contain one derivative of V, and therefore 



cannot be bounded in the norm. However, for N' < N, we can obtain an estimate for the norm by 



v 



combining Proposition A. 10 (7.3.271 and the uniform bound provide by the constant K. The inequality we 
obtain is 



(7.3.28) 



\ H „. iT , <C||V(0)||^y' /Ar .exp(CT'), 



where in (7.3.28) , C = C(N',N,0 2 ,K). 

Remark 7.3.2. The estimate ( 7. 3. 27) is a limiting version of the "conical" estimate 
(7-3.29) l|V(t)|| H *-i (Et(r _ t) < C£||V(0)||**- 1(EOir) •exp(Ct), 



where we are using notation defined in Lemma 7.2.2 A proof of ( 7.3.29 1 can be constructed using arguments 



similar to the ones used in our proof of (7.2.301. Inequality (7.3.291 shows that two solutions that agree 
on Eo,r also agree on £ tjT ._ t . By translating the cone from Lemma |7.2.2 so that its lower base is centered 
at the spacetime point x, we may produce a translated version of the inequality. Thus, we observe that a 
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domain of dependence for x £ M. is given by the solid backward light cone in M. with vertex at x; i.e., the 
past (relative to x) behavior of a solution to the EOV outside of this cone does not influence behavior of 
the solution at x. Similarly, a domain of influence of x is the solid forward light cone with vertex at x; the 
behavior of a solution at x does not influence the future (relative to x) behavior of the solution outside of 
this cone, a fact which justifies our claim made in Section |5.4| that the fastest speed of propagation in the 
EN K system is the speed of light . In [6 , Christodoulou gives an advanced discussion of these and related 
topics for hyperbolic PDEs derivable from a Lagrangian. 



This completes our abbreviated proof of Theorem [T] 



□ 



7.4. Proof of Theorem [2j We now provide a detailed proof of Theorem [2j 

7.4.1. The setup. Let {V m } b e the sequ ence o f initial data from the hypotheses of Theorem bl converging 



in H$ to V 



7.0.1 



and 



By corollaries 
solutions V™ and V respectively to (4.1.1 
the property 



7.0.2 



for all large m, the initial data V™ and V launch unique 



(4.1.261 that exist on a common interval [0, T'] and that have 



V([0,T'] xK 3 ) C 2 ,V 
K = K(N,0 2 ,\\WV\\h. 



l ([0,T'] x R 3 ) C 02- Furthermore, for all large to, with 



(7.4.1) 



V 



\+i,A, S), we have the uniform (in m) bounds 
Hirw-i xi HI V™ || \ H N T t, HI d t V r 



Iff™ T' 



Iff' 



< K, 



where [0, T] is the interval of existence for V furnished by Theorem [T] In this section, we will show that for 
all large to, V™ exists on [0, T] and that 

(7.4.2) lim HI V m - V \\\ h n T = 0. 

m — >oo 

The proof we give here is inspired by a similar proof given by Kato in |18j . We use results and terminology 
from the theory of abstract evolution equations in Banach spaces, an approach that streamlines the argument. 
We also freely use results from the theory of integration in Banach spaces; a detailed discussion of this theory 
may be found in [32]. We begin by rewriting the linearization of the EN K system around V and V m as 
abstract evolution equations in the affine Banach space H^(R 3 ). In this form, the linearized systems are 
respectively written as 



(7.4.3) 
(7.4.4) 



d t Z + A(V)Z = f (V) 

a t z + ^(v m )z = f(v m ), 



where f(-) is a smooth function on O, and f(V) = 0. Here, the symbol Z stands for all 10 components of 
a solution to a linearized system, and the operator A(-) is a first order spatial differential operator with 
coefficients that depend smoothly on its arguments. We state for clarity that the first 5 components of the 
inhomogcncous terms f(V m ) are given by (^°(V m )) _1 • (S(V m ), 6(V m ), • • • 7 ^ (V m )) Trans P ose , where the 
matrix- valued function A°(-) is defined in ( |4.2.20 1 and the scalar- valued functions J, ©, • • • $) 3 arc defined in 
(14.2. Ill) - (14.2.131). 



We will make use of the pseudodiffcrcntial operator 

dcf/-, (1 JV/2 



(7.4.5) 



S = (1 - A) 1 



which is an isomorphism between H and L ; i.e., S S C(H , L ) and 
S- 1 G C(L 2 .H N ). 

7.4.2. Technical estimates. In this section, we provide some technical lemmas that will be needed in our 
proof of Theorem [2] For certain function spaces X, there exist evolution operators 

(7.4.6) U(t, i'), U m (t, t'):X^X 

defined on A T , = {Q <if <t< T'} that map solutions (belonging to the space X) of the corresponding 
homogeneous version of the linearized systems (7.4.31 and ( |7.4.4 | at time t 1 to solutions at time t. The 
relevant spaces in our discussion are X — L 2 and X = H N . In the following three lemmas, we describe the 
properties of the operators U(t,t') and U m (t,t'). Complete proofs are given in [TB], [17], and [IB]; rather 
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than repeating them, we instead attempt to provide some insight into how the proofs relate to the methods 
described in this paper. 

Lemma 7.4.1. U(-,-) and U m (-,-) (for m > 0) are strongly- continuous maps from Ay into C(L 2 ) n 
C(H N ). Furthermore, there exists a C(K) > such that ||| U a n\\\ M m IIIl 2 a M \\\h n a , 
,\\\U m \\\ HN ,,<C{K). 



Remark 7.4.1. Lemma 7.4.1 is essentially a consequence of the fact that the uniform bound (7.4.1 1 for V 

and V m allows for uniform Sobolev estimates to be made on the L 2 (or H N ) norm of L 2 (or H N ) solutions 
to the linearized equations. 

Remark 7.4.2. By the regularity result of V furnished by Theorem [l] Corollary |A.7j Remark A. 5 and 
06X2) , the right-hand side ( [7X3| is an element of C°([0, T],H N ). Given a function Z € it follows 
from Lemma 7.4.1 and standard linear theory (via Duhamel's principle) that there exists a unique solution 
Z E C°([0, T], Hy) to (7.4.31 with initial data equal to Z. An analogous result holds for solutions to (7.4.41. 

Lemma 7.4.2. I4 m (t,t') converges to U{t,t') strongly in C(L 2 ) as m — > oo. Furthermore, the strong con- 
vergence is uniform on Aj-i. 

Remark 7.4.3. By smoothing the initial data, a solution Z E C°([0, T],L 2 ) to either d t Z + A(V)Z = 
or d t Z + A(V m )Z = f(V m ) can be realized as the limit (in the norm ||| • |||x,a y) of a sequence {Z fe } C 
C°([0, T], H N ). Therefore, to prove Lemma 



7.4.2 



one only needs to check that given initial data Z e H 



N 



we have that 
(7.4.7) 



lim 

m — >oo 



(U m (;0)-U(;0))Z 



\l 2 ,t- 



0. 



Based on Lemma 7.4. 1| and (7.3.281, which can be used to show that for N' < N we have V r 
C°([0,T'],H^') n Cl([0,T'] xl 3 ), ( 7X7] follows from the method of energy currents. 



V in 



Lemma 7.4.3. There exist operator-valued functions B : [0, T] — > C(L 2 ) 
B m : [0,T'] -> £(L 2 ) smc/i ttat 

S^(V(i))S- 1 = >t(V(t)) + B(t) 
S^(V m (t))S" 1 = ^(V m (<)) + 6 m (t). 



(7.4.8) 
(7.4.9) 

Furthermore, for all t' with < t' < T' , B and B m satisfy the estimates 
(7.4.10) HI B m - B \\\ L 2 tt , < C(K) HI V m - V 



\H N ,t' 



(7.4.11) 



and HI B 



B n 



lL2jT , < G{K). 



Remark 7.4.4. A modern treatment of commutator estimates that can be used to prove Lemma |7.4.3| is 
located in [33 . 

Our proof of Theorem [2] also requires the following lemma, whose simple proof is based on Duhamel's 
principle: 



Lemma 7.4.4. Let Z E H^., and let Z E C°([0, T], H^) denote the unique solution to (7.4.31 with initial 

data equal to Z as furnished by Remark 7.4-2 Then for 
t E [0,T], S(Z — V) satisfies the Duhamel formula 



(7.4.12) 



S(Z(*)-V)=W(t,0)S(Z-V)- / U(t,t , )B(t')S(Z(t , )-V)dt' 



U{t,t')Sf(V{t'))dt'. 



An analogous result holds for the linearization of the EN K system around V 1 



Proof. We apply S to each side of the equation satisfied by Z — V and use Lemma 7.4.3 to arrive at the 
equation 

(7.4.13) d t [S(Z - V)] + A(V)S{Z - V) = -B(V)S(Z - V) + Sf (V). 
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Thus, S(Z — V) is a solution to the same linear equation that Z solves, except the inhomogeneous terms for 
S(Z — V) are given by the right-hand side of (7.4.131 and the initial data are given by S(Z — V). Equation 
(7.4.12) now follows from D uham el's principl e. N ote that Sf(V) is well-defined since f(V) = 0, and we can 
therefore apply Proposition |A.8 and Remark A. 6 to conclude that f(V) S H N . □ 

7.4.3. Proof of Theorem^ We will now demonstrate ( |7.4.2 1 by providing a proof of the equivalent statement 
(7.4.14) lim HI S(V m - V) ||| L 2 T =0. 

m — >oo ' 

Lemma |7.4.4 implies the following equality, valid for < t < T' : 

S(V ro (t) - V(t)) = U m {t, 0)S(V m - V) + (U m (t, 0) - U{t, 0))S(V - V) 

U(t, t')B(t')S(V(t') - V) dt' - f U m (t, t')B m (t')S(V m (t') - V) dt' 



(7.4.15) + / U m (t,t')Sf(V m (t'))dt' - / U(t,t')S{(V(t'))dt'. 

Jo Jo 

By Lemma [7.4.1 1 we have that 

(7.4.16) HI U m (t, 0)S (V m -V) \\\l*,t><C(S,K) \\ V m -V|| H » . 
Furthermore, if we define 
(7.4.17) 



sup \\(U m (t 7 0)-U(t,0))S(V~V)\\ L 2 

te[o,T'] 



then Lemma |7. 4. 2| implies that 
(7.4.18) 



lim dr, 

m — ^oo 



We now rewrite the second line of (7.4.151 as 



(U(t, t') - U m (t, i'))B(*')S(V(i') - V) dt 1 



n 



+ / U m {t,t'){B{t')-B m (t'))S{V{t')-V)dt' 



(7.4.19) 



+ / U m (t,t')B m (t')S(V(t')-V m (t')) dt'. 



By (7.4.1 1, Lemma 



7.4.1 



and Lemma 



7.4.3 



for < t < Tl < T' the L 2 norms of the second and third 



integrals in (|7.4.19|) are each bounded from above by 
C(S, K)T» 



"V m — "V | | \ f{N y 

We simila rly sp lit the third line of (7.4.151 into two terms and use ( 7.4.1} , Lemma 7.4.1| Proposition A. 8 



A.6 



to bound the L 2 norm of one of them, namely /„' U m {t, t')S(f(V m (t') - f(V(t'))) dt', from 



and Remark 
above by 
C(S,K)T* HI V m -V\\\ H »,T* ■ 

Combining these estimates with (7.4.161 and (7.4.171, we take the L 2 norm of each side of (7.4.151 followed 
by the sup over t £ [0,T*] to arrive at the inequality 

HI S(V m -V) IHi, , T . < C(S,K)\\V m ~V\\ HN +d m 
+ C{S,K)T* HI V m -V \\\ h n iT , 

sup \\{U{t, t') - U m (t, t'))B(t')S(V(t') - V)|| L 2 dt' 
o te[o,r.] 

(7.4.20) +/* sup \\(U m (t 7 t')-U(t,t'))Si(V(t'))\\ L 2dt'. 

Jo t£[0,T»] 

We now choose small enough so that 



(7.4.21) 



C(S,K)T* HI v m -V|||^,r.< 2 ll|S(V m -V) m^r. 
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from which it follows that 

HI S(V m - V) 



(7.4.22) 



| L » iT . <2C(S,if)||V m - V\\ HN +2dm 

sup || {U m {t,t') -U{t,t'))B{t')S{V{t') -V)\\ L * d£ 
11 te[o,T„] 

sup ||(W m (t,0-^(i,i'))S(f(V(i'))-f(V))|| i2 di', 

te[0,T„] 



where in (7.4.221, we have used the fact that f(V) = 0. 

Lemma 7.4.3| Proposition A. 6 Remark A. 5 Proposition A.8[ and Remark 



By (7.4.1 1, Lemma 7.4.1 



A. 6 the two integrands in (7.4.221, viewed as functions of t' , are uniformly bounded by C(K) on [0,T*]. 
Furthermore, by Lemma 7.4.2 the integrands converge to pointwise in t' as to — * oo. Therefore, by the 
dominated convergence theorem, the two integrals in (7.4.221 converges to as to — > oo. Recalling that by 

V m - V|| 



hypothesis we have that lim,, 
(7.4.23) 



lim 



\gn = 0, and also using (7.4.181, we conclude that 

S(V m -V) ih^tv 



0. 



To extend this argument to the interval [0, 2T*], let e > and choose Too large enough so that m > too 
implies that ||| V m — V |||hw 1 t«< e /(4C(S, K )) . Starting from time T*, we may argue as above to show that 



(7.4.24) 



lim sup sup 

ro— >oo te[T„,2T» 



|S(V m -V)||^ 



1 



Thus, we can choose toi > too such that ||| S (V m —V) |||z, 2 ,2T,< c when to > m\. Continuing in this 
manner, we may inductively extend this argument to the interval [0, T'\. We state for emphasis that the size 



of T* required to satisfy the inequality (7.4.211 depends only on C(S,K). Consequently, the length of the 
time interval of extension may be chosen to be the same at each step in the induction. 

We now show that this argument can be extended to the entire interval [0,T] on which V exists. Define 



T m ax = sup{T' | V and the V m exist on the interval 
[0, T'] for all large to and lim 1 1 1 V 1 



- V 



\H N ,T 



= 0}. 



We will show that the assumpti on T rn ax < T leads to a contradiction. 

By Theorem [l] and Corollary 7.0.1 for each t £ [0,T], there exist an H N neighborhood Bg t (V(t)) of V(i) 
with positive radius St and a A t > such that initial data belonging to B$ t (V(t)) launch a unique solutiorpj 
that exists on the interval [t, t + A t ] (the term "initial" here refers to the time t). By continuity, V([0, T]) is a 
compact subset of . Therefore, there exist <5 > and A > such that initial data belonging to Bg(V(t)) 



launch a unique solution that exists on the interval [t, t+ A]. Furthermore, by Corollary 7.0.2 there exists a 
C un iform > such that for any "initial" data V contained in the ball _Ba(V(i)), the corresponding solution 
V to the EN K system satisfies the bounds 



V 



<a 



uni f orm 



(7.4.25) 



\{t,t+A] 



<c, 



uni f orm • 



We emphasize that S and A are independent oft belonging to [0, T], and that C uniform is independent of 
the data. Note that as a consequence of this reasoning, it follows that V exists on the interval [0, T + A]. 

The contradiction is now easily obtained: assume that T max < T. Then according to the above paragraph, 
initial data belonging to Bs(V(T max — | A)) launch a solution that exists on the interval [T ma x — \ A, T max + 
~A]. Furthermore, for all large to, ~V m (T max — |A) is contained in Bs(V(T max — |A)). Therefore, for all 
large to, V m can be extended to a solution that exists on [0, T max + | A] . In addition, using ( 7.4.25 ) , we have 
that for all large to, ||| V m \\\ h n TmQx + i A < G uni f orm , and ||| <9 t V m Hl^w-i r m „+iA< C U niform- Therefore, 

V 2 ' 2 

we can repeat the entire argument given in Section 7.4 substituting T max + ^A in (7.4.1 ) in place of T' and 
T, and C un iform in place of K, to show that linim^oo ||| V™ — V |||^n r mal +iA : 
definition of T max and completes the proof of Theorem [2j □ 



0. This contradicts the 



m This solution may escape 02, but this is not a difficulty since O2 <E O; the solution still has some "room" 
remain in a compactly supported convex subset of O . 



left to evolve and 
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A. Appendix 

In this Appendix, we use notation that is as consistent as possible with our use of notation in the body of 
the paper. To conserve space, we refer the reader to the literature instead of providing proofs: propositions 



A. 6 and A. 8 are similar to propositions proved in chapter 6 of |15j . while Proposition A. 9 is proved in [20 



The corollaries and remarks below are straightforward extensions of the propositions. With the exception of 
Proposition |A.10| which is a standard Sobolev interpolation inequality, the proofs of the propositions given 
in the literature are commonly based on the following version of the Gagliardo-Nirenberg inequality [27 , 
together with repeated use of Holder's inequality and/or Sobolev embedding: 

Lemma A. 5. If i, k G N with < i < k, and V is a scalar-valued or array-valued function on R d satisfying 

V G L°°(R d ) and || V( fe )V|| L 2 (Rci) < oo, then 

(Al) ||vWV|| L2fcA < C(£:)||V||^||V( fe )V||! 2 . 

Proposition A. 6. Let O2 C K" be a bounded open set, and let j, d G N with j > |. Let V : R d — > E" be an 

element o/-£P(R d ), and assume that V(R d ) C 02- Let F G C^iO-i) be a q x q matrix-valued function, and 
let G G Hi (R d ) be a q x q (q x 1) matrix-valued (array -valued) function. Then the q x q (q x 1) matrix-valued 
(array-valued) function (FoV)G is an element of iP(R d ) and 

(A.2) \\(F°V)G\\ Him < C(j,d)\F\ j!02 (l + \\y\\ j H , {Md) )\\G\\ Him . 

Corollary A. 7. Assume the hypotheses of Proposition \A~d\ with the following changes: 

V, G G C°([0, T], _ff J (R d )). Then the q x q (q x 1) matrix-valued (array-valued) function (F o V)G is an 
element of C a ([0,T],W(R d )). 

Remark A. 5. We often make use of a slight modification of Proposition |A.6| in which the assumption 

V G i/ J (R d ) is replaced with the assumption V G H^.( W. d ), where V G R" is a constant array. Under this 
modified assumption, the conclusion of Proposition |A.6| is modified as follows: 

(A.3) \\{FoV)G\\ Hi <C(j,d)\F\ jt62 (l + \\V\\^,)\\G\\ Hi . 



A similar modification can be made to Corollary |A.7[ 

Proposition A. 8. Let O2 C R" be a bounded open set with convex closure O2, and let j, d € N with j > | 
Let F G C\ (O2) be a scalar or array-valued function. Let V, V : R d — > R™, and assume that V, V G ZP (1R ) 
Assume further that V G 2 and V(R d ), V(R rf ) C 2 . Then foV-foVe W(R d ) and 

(A.4) \\FoV~FoV\\ H3 <C(i,d,||V||^,||V|U 3 )|^ +1)a J|V-V|| H i. 

we may replace the hypotheses V, V G H : >(W. d ) from Proposition 



A.5 



A.8 



with 



Remark A. 6. As in Remark 
the hypotheses V, V G (R d ) , where V is a constant array, in which case the conclusion of the proposition 
is: 

(A.5) \\(FoV)-(FoV)\\ Hi <C(7,d,||V||^,||V||^)|F| i+li< 5 a ||V-V|| ff ,. 

Furthermore, a careful analysis of the special case V = V gives the bound 

(A.6) \\F oV-FoV\\ w <C(j, d)\dF/dV\ ] _ h02 (l + ||V||J£)(||V||^), 

in which we require less regularity of F than we do in the general case. 

Proposition A. 9. Assume the hypotheses of Proposition \A~(j\ with the following two changes: 
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(1) Assume j > | + 1. 

(2) Assume that G G ff^QR"*). 

Lei fc G N 1 < fc < j, and let a be a spatial derivative multi-index with \a\ = k. Then 

\\d 3 ((FoV)G)-(FoV)d s G\\ L 2 
(A.7) < C&QldF/dVl^&iWVWn + ||V||k)||G|| H ,-i. 



we may replace the assumption V G H 3 (R ) in Proposition A. 9 with the 



Remark A.7. As in Remark A. 5 

assumption V G H^.(M. d ), where V is a constant array, in which case we obtain 

\\d s ((F oV)G) - (F oV)d s G\\ L 2 
(A.8) < C&QldF/dVl^o^WVW^ + ||V||^)||G|| ff i-i. 

Proposition A.10. Let N',N eR be such that < N' < N, and assume that F G H N (R d ). Then 
(A.9) \\F\\ HN ,<C(N\d)\\F\\l 2 N ' /N \\FCT. 
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